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ABS TRACT 

In this work the coaxial waveguide antenna is treated by the 
Wiener-Hopf technique and the transient radiation from a cylindrical 
monopole is developed in the light of the rigorous results obtained 
from the Wiener-Hopf analysis. Analytic expressions are derived for 

(1) the electromagnetic fields in the feed line and (2) the far zone 
radiation fields of the coaxial waveguide antenna, with time harmonic 
excitation voltage. Complete characterization of the transient beha- 
vior is also found for (1) the fields interior to the feed line and 

(2) the radiated fields for excitation voltages arbitrary in their 
time dependence to the extent that kb,ka « 1 . This corresponds 
to the case of a thin antenna and excitation voltage with a non-zero 
rise time, specifically chosen so that frequencies violating the 
restriction kb,ka « 1 are negligible. 

The transient radiation from the cylindrical monopole is 
developed in a closed analytic form which is relatively easy to 
interpret and apply. The expressions found offer an alternative to 
transient analysis by conventional methods requiring numerical tech- 
niques involving extensive computer calculations. They are also the 
basis for an uncomplicated procedure to synthesize a desired behavior 
of the transient radiation from cylindrical monopole antennas. 
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1. INTRODUCTION 

This dissertation presents an analysis of the coaxial waveguide 
antenna by the Wiener-Hopf technique and the determination of transient 
radiation from the cylindrical monopole antenna. This work is the 
first successful treatment of a semi-infinite coaxial structure by the 
method of Wiener-Hopf and the first rigorous determination of transient 
radiation emanating from an open ended semi -infinite waveguiding struc- 
ture. In addition, a new and simple method for calculating the 
transient electromagnetic radiation of a cylindrical monopole antenna 
is presented. 

The study of electromagnetic boundary- value problems by the 
Wiener-Hopf technique began in the early 1940's. Among the most notable 
of the early workers were Schwinger [1] and Copson [2]. The early 
applications of this powerful method were inthe rigorous solution for 
the fields within a bifurcated parallel plate waveguide and diffraction 
of plane waves by a conducting half-plane. The method was later used 
to solve the problem of radiation impinging on or emerging from the 
open end of an unflanged circular pipe and a pair of parallel plates. 
Much of the early history of application of the Weiner-Hopf method is 
detailed by Bouwkamp [3] and an authoritative account of the method was 
published in 1958 by Noble [4]. In this work the integral equation 
approach of the early workers is not stressed . A more direct simple 
method due to Jones [5] is emphasized. It is the application of 
Wiener-Hopf technique due to Jones that will be employed within this 
report. 
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There have been several recent investigations of the transient 
behavior of antennas [6-15]. The focus of many of these investigations 
is on the time behavior of the radiated electromagnetic field. Refer- 
ences [6,8] and [14] present theoretical as well as experimental 
results. In most of the prior research efforts, the problem is formu- 
lated in terms of circuit theory and the antenna is treated as a lumped 
circuit element. The feedline and the exciting sources are replaced by 
their Thevenin or Norton equivalents. The formulation of the problem 
in terms of circuit theory yields complex frequency domain fields which 
must be Fourier transformed if the time behavior is of interest. Such 
calculations require numerical methods due to their extreme complexity. 

The work of Schmitt, Harrison and Williams [6] is a determina- 
tion of the transient radiation and reception performance of a thin, 
finite length cylindrical monopole over a perfectly conducting infinite 
ground plane with a coaxial feed. In the circuit analysis formulation 
(frequency domain) the antenna was represented by an impedance and the 
radiated field related to the voltage across the antenna impedance 
through a proportionality factor, called effective height. The time 
dependent behavior is obtained by Fourier transforming the frequency 
domain field. A similar procedure was used to obtain the properties of 
the antenna in reception; the induced voltage on the receiving antenna 
was related to the incident electric field intensity through the effec- 
tive height. The theoretical work of Schmitt, Harrison and Williams 
was supplemented in 1971 by Abo-Zena and Beam [9] to include the time 
behavior of near zone field intensities and in 1970 by Palciauskas and 
Beam [10] to include the far zone time behavior at latitudinal angles 
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other than 90°, which was the only angle considered by Schmitt et al . 

King and Harrison [8] formulated a frequency averaged reflec- 
tion coefficient to theoretically obtain the reflected electromagnetic 
fields on the coaxial feed line of an infinitely long cylindrical 
antenna or for the time interval before the excitation pulse can 
reach the tip of a finite length antenna. Use of frequency averaging 
is necessitated by the fact that the antenna input admittance is 
frequency dependent, which constrains the conventional transmission 
line reflection coefficient to be frequency dependent. 

The transient current on a circular tubular infinitely long 
antenna that is excited by a voltage which is a step function of 
time was determined in 1961 by T. Wu. He later (1969) furnished a 
corrected analysis for this current in Reference [11]. An earlier 
work (1960) by Brundell [12] treated the more general problem of 
determining the space and time behavior of the electromagnetic field 
for such an antenna. This work was also supplemented in 1970 by 
Latham and Lee [13] who found the early and late time behavior of 
the field radiation from a hollow infinite cylindrical antenna 
excited by a voltage that is a step function of time. 

Lamensdorf [14] determined that the time behavior of the 
coaxial cone antenna in reception is primarily a time derivative of 
the incident electric field intensity falling on the aperture. This 
work also contains a heuristic discussion of the time dependent elec- 
tric field radiated by a cone antenna. Chang [15] analyzed the tran- 
sient reception characteristics of the annular slot antenna using the 
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circuit theory approach. References [14] and [15] are recent works, 
1970 and 1971 respecti vely, and have application as part of a receiv- 
ing system for measuring scattered impulse responses of various 
objects or as a receiving antenna on an aircraft or guided missile. 

In t he systematic study of radiation from the coaxial aperture 
antenna which follows, analytic expressions for the electromagnetic 
fields of the coaxial aperture antenna are given. These expressions 
account for waves traveling in both the forward and reverse directions 
on the feedline and the radiation fields for harmonic time dependent 
excitation. Complete characterization of the transient behavior is 
found for the fields interior to the feedline and the radiated fields, 
for excitation voltages which are arbitrary to the extent that 
kb,ka « 1 . This corresponds to the case of a thin antenna and 
excitation voltages with non-zero rise times, specifically chosen so 
that frequencies which violate the restriction kb,ka « 1 are neg- 
ligible. 

The annular slot antenna is considered to be a special case of 
a cylindrical monopole antenna (h=0). The fields radiated from the 
annular slot antenna are approximated by well known methods in which 
the incident field on the aperture is considered to be the aperture 
field. The radiated fields so determined are interpreted in light of 
the rigorous results obtained from the coaxial waveguide antenna and 
a model is obtained for a cylindrical monopole antenna in which the 
height of the radiating element is greater than zero (h >0). The 
transient radiation field intensities are obtained from this antenna 
model. The determination of the transient radiation through the 
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procedures developed in this research offer an alternative to comput- 
ing transient behavior by numerical methods where extensive machine 
calculations are required. The procedures reported on herein can also 
be used as the basis of an uncomplicated synthesis procedure to produce 
a desired time behavior in fields radiated from the coaxial waveguide 
and the cylindrical monopole antennas. 

The results of this research should prove to be very useful in 
the areas of 1) determination of antenna properties (driving point 
impedance and radiation pattern) through time domain measurements [16], 
2) identification of and discrimination between conducting bodies 
through radiation scattering from their surfaces and 3) the calibration 
of equipment designed to measure the effectiveness of shields designed 
to protect equipment or installations from intense, short rise time 
electromagnetic pulses. In each of the above areas a precise knowledge 
of the time-behavior of the radiated electromagnetic field is required. 
Additionally the ability to synthesize a particular time behavior in 
the radiated field is highly desirable and very important. 
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2. THE COAXIAL WAVEGUIDE: AN APPLICATION OF THE WIENER-HOPF TECHNIQUE 

In this chapter the frequency (spatial) domain electromagnetic in- 
tensities for an open ended coaxial waveguide will be derived. This will 
be done using the methods of Wiener and Hopf [4]. Use of this technique 
will provide a determination of transverse electromagnetic mode (TEM) 
as well as higher order mode intensities within the structure. 

We take the coaxial waveguide structure to be that of Figure 2-1. 
The inner conductor is a hollow tube. Both outer and inner conductors 
are assumed to be semi -in finite in length, infinitesimally thin and 
perfectly conducting. The cross section is circular providing for a 
coaxial annular region between the conductors. We also assume that the 
waveguide is embedded in a homogeneous, isotropic, dielectric medium 
having the constitutive parameters of free space. 

Specification of the Problem 

We begin the specification of the problem by listing the Maxwell 
equations: 

VXH = J + lf 

V • D = p 

V • B = 0 (2.1) 

The fields which satisfy the above equations must also satisfy the con- 



stitutive relations: 
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Figure E-1. The coaxial waveguide antenna 
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where 




Under the assumption that the coaxial waveguide has been excited in 
TEM mode, there is a traveling electromagnetic disturbance proceeding 
in the positive z direction. The time origin is taken to be the 
point in time when the leading edge of the TEM mode exciting field is 
incident upon the aperture located in the z = 0 plane. 

The specific method used to initially excite the coaxial struc- 
ture and the field distribution prior to the time origin are not of 
interest in this work. 

It is well known that the TEM mode field of a coaxial waveguide 
is independent of angular variation [17, p.326] in the cylindrical 
coordinate system. The coaxial structure also does not introduce any 
boundary conditions which have angular dependence. From these observa- 
tions we conclude that the scattered field must also be characterized 
by the same symmetry (i.e., independent of angular variation in cylin- 
drical coordinate system). 

In cylindrical coordinates the field components can be written 



as 



E 




and 



H 




where a ,a.,a is the triad of unit vectors along the coordinate axes, 
p’ (j)’ z 
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We now Fourier transform the Maxwell equations to change from a 
time domain to a frequency domain representation. The Fourier trans- 
form pair [18] is defined as follows: 



F(w) = — f f(t) e^“^ dt 



f(t) = ^ 



F(w) do) 



When the first two Maxwell equations are multiplied by 
operated on by the integral operator, the results are: 

oo oo 







V X E e 



iwt 



dt 




f 3B ^iwt 

I H 



dt 



1 giwt 
^/2^\ 



and 



and 



1 



I y X 



H e'^^t 



dt = 



1 






dt 



Interchanging the curl operator and the integral sign on the left hand 
side and integrating by parts on the right hand side yields equations 
in the frequency domain. To form these results, use has been made of 
the fact that the field intensities must be zero as |t| ->■<». 

V X E(w) = + iw B(w) 

y X H(w) = J(w) - iu) D(o)) 



It similarly follows that 







■■■ ■ ■ ■ ■ 







^ I 
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V • D(w) = p(w) 

V • B(w) = 0 

B(w) = li^H(a)) and D(w) = e^E(o)) 



When the field components in cylindrical coordinates are substituted into 
the frequency domain Maxwell equations for a source free region (i.e., 
J(w) is set to zero), we obtain the following relations between the 
field components; 



iw B 

P 

iw B^ 
iw B, 



_ —± 
dz 

8E, 9E, 

_1 + -1 
P 9p 

9E 9E, 

_£. I 

9z 9p 



( 2 . 1 . 1 ) 

( 2 . 1 . 2 ) 

(2.1.3) 



and 





(2.1.4) 


P 9p 


(2.1.5) 


9H 9H, 

1 

9Z 9p 


(2.1.6) 



By using the constitutive relations and combining equations (2.1.3), 
(2.1 .4) and (2. 1 .5), we derive the scalar wave equation for , 



iw H. 
o <p 



‘ Iz li ”4.^ • If 



This equation* upon simplification, becomes: 
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H^(.) 



0 



( 2 . 2 ) 



where ye = -o 
o o 

c 

produce 




IL. 

P 9P 



Similarly equations 





2 

0 ) 1 
T-l 

c 



E^(.) = 



( 2 . 1 . 1 ). 



0 



(2.1.2) and (2.1.6) 



(2.3) 



From the foregoing equations it is apparent that the solution of 
(2.2) and (2.3) for and is enough to completely determine all field 

components. The work incorporated in Appendix A shows that E^ must 
be zero for the coaxial waveguide, therefore the only field intensity 
components present are E^ and E^ • The rest are identically 
zero. The analysis has been reduced to solving the boundary value 
problem 

2 2 

r ^ H,(p,Z,0)) = 0 (2.4) 

9Z^ 9p^ P ‘P 



H^(0,z,a)) = 0 ; -« < z < « (2.4.1) 

E^(a,z,to) = E^(b,z,to) = 0 ; Z£0 (2.4.2) 



The Sonmerfeld radi ati on condi ti on , 

3H. 

lim r(:^ - ik H.) = 0 must also be satisfied. 



At this point it is useful to point out that in the above equations 




(2.4.3) 



i ^ 
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is a scalar quantity and equation (2.1.5) can be written as 




(2.5) 



This will be of use later. 

The well known expressions for the TEM mode fields for the 
coaxial waveguide [17] will be given here to set the notation. We 
take the field intensities incident on the aperture from the negative 
z direction to be 



v(w) is the frequency domain voltage between the inner and outer con- 
ductors . 

The space containing the waveguide is divided into three regions 
portrayed in Figure 2-2. Equation (2.4) applies to each of these regions 
The field satisfying this equation in region one is designated 



H (p,z,co) for region three. The equation for each of the regions is 
a partial differential equation in the variables p and z . By 
Fourier transforming in the spatial dimension z , an ordinary differen 
tial equation in the variable p is obtained. 

Following the notation of [4] , the spatial Fourier transform 



, vfg.) 



( 2 . 6 ) 



^ lx\ b/a ^ 




(2.7) 




pair is defined as follows: 



-13- 



Region I 




Region 1 b £ p 

Region 2 a £ p £ b 

Region 3 p £ a 



Figure 2-2. Division of the space surrounding the coaxial waveguide 
antenna into regions 



F(a) 
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f(z) dz 
F(a) e da 



The "plus" and "minus" functions are defined as 



F''(a) = ^ 



f{z) e^“^ dz 



0 



F"(a) = -L 



✓57 J 



f(z) e^“^ dz 



where 



F(a) = F'^(a) + F‘(a) 



( 2 . 8 ) 

( 2 . 8 . 1 ) 



( 2 . 8 . 2 ) 



(2.8.3) 



(2.8.4) 



a = a + ix is complex . 



Using well known theorems involving Fourier transforms in the 
complex plane [4], the regions of analyticity of F^(a) and F"(a) 
will be found. Suppose that we are given x_ < r < t_j_ , all real 
niinbers. In addition suppose that 



|f(z)| < A exp(x_z) as z ->■ 



where A is a constant. Then 

M 



F'"(a)| < 

✓57 



f(z) e'*^" dz 



0 



A 



(ia + x_)z 



dz 



where M < 



oo 



is chosen large enough for the inequality to be valid. 
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The last integral of the above expression converges if and 

only if 

(ia+ T )z 

e ~->0 asz->-“ 

or 

(io - T + T_ )z 

e ->-0 as z 

This last expression is valid provided 

T > T_ 

By similar reasoning, it can be shown that if |^"(z)| < B exp(x^z) 
as z -»• -00 (B constant) then |F”(a)| is convergent provided 

T < 

Therefore F^(a) and F'(a) are convergent in the strip 
T_ < T < of the a plane. 

We also note that if |F^(a)| and |F~(a)| are convergent 
in a region of the complex plane, then the magnitude of nth ordered 
derivative (n < oo) of F^(a) and F"(a) with respect to a also 
exists in that region. In fact such functions are nth order continu 
ously differentiable in the region. 



^ (F^c<)) 



da' 






da" 




^n 



M 






z"f(z) e'“" dz 



OO 




(ia + T_)z 
e 



< 



dz 
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where M < ~ is properly chosen to ensure the inequality. Concen- 
trating on the term of significance, we write 



z" e 



(ia+T)z (ia+T)z 

■ dz = z" e 



M 



1 a + T 



, (ia+T )z 
z""' e ■ dz 



After n integrations by parts, the above becomes 

1 «i (ia+T )z “ „ „/ , (ia + T )z 

V » (z) (-1) p - + { i^n n(n-l)--»l ^ 



P=1 ' (ia+T ) 



(ia+T )' 



which is finite for finite n provided 



T > T_ 

Again by a similar discussion, it can be shown that F”(a) is nth 
ordered continuously differentiable in the open region t < t^ provided 
|f(z)| < Be^+^ as z • 

Since all that is required for analyticity in a region of the 
complex plane is that the first derivative of the function exist in that 
region, we are assured that F_j_(a) and F"(a) are analytic in the 
strip T < T < T^ of the a plane provided of course that 

T_Z 

I f(z) I < A e ” z “ 

and 

I f(z) I < B e z 



Upon multiplication of equation 2.4 by e’^^ 

c“ yzi\ 

dz , the following equation in 



ing with the integral 



and operat- 



p emerges. 



—OO 
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P 9p 



- -^ + k^] H ,(p,a,w) = 0 
P ^ 



Integrati ng by parts twice and imposing the conditions that H^(p,z,co) 
and ■^H^(p,z,to) -^0 as jz| , produces an ordinary differential 
equation in p , 




(2.8.5) 




( 2 . 8 . 6 ) 



It is noted at this point that the incident magnetic intensity does 
not satisfy the differential equation (2.8.5) since it does not converge 
in the limit as z -^ -«> . This will be explored in detail a little 
later. 

For region 1: 




(2.9) 




as p -> «> 



(2.9.1) 



^ [p H-^(p,a,io)] 



p=b 



0 



(2.9.2) 



Equation (2.9.1) follows from the Sommerfeld radiation condition and 
equation (2.9.2) follows from equation (2.5) and the boundary condition 
(2.4.2). 



-18- 



For region 3: 



[ 



+ 

^ p *. 



1 2 

“Y ■ Y ] H , (p,a,io) 
P 



( 2 . 10 ) 



H (0,a,io) = 0 

$3 



( 2 . 10 . 1 ) 



^ [p H'^(p.a.(o)] 



= 0 



( 2 . 10 . 2 ) 



p=a 



Equation (2.10.1) and (2.10.2) result from imposing (2.4.1) and (2.4.2) 



For region 2: 

In this region the magnetic intensity will be written as a 

linear combination of the incident TEM magnetic intensity and another 

term ifj (p.z.io) . 

$2 



H^^(p,z,io) + 



ip (p.z.io) 
^2 



( 2 . 11 ) 



The term tp (p,z,to) represents obviously all of the reflected fields that 
^2 

exist within the feedline, on z < 0 . For the other extreme of 



z > 0 , ip. (p,z,w) must account for the fields which exist in free space 
4>2 

from sources within the feedline. It must be representable as 



ip^ (p,z,w) = (-h 1 + other field terms) , z > 0 
4>2 4 > 

since cannot be present in the half space z > 0 because it vio- 
lates requi red boundary conditions. 
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Since the intensity H, (p,z,co) must be a solution to the dif- 

^2 . 

ferential equation (2.4) the terms and (p,z,co) must satisfy 



the following; 



and 



[- 



3Z 

.2 



9p 

.2 



_+ 

2 p 3p 



V + k ] = 0 

. ^ Q) 



r 9 , 9 , 1 9 



^ + k^] \p (p,z,o)) 

p ^2 



= 0 



The incident field identically satisfies the 

first of the above equations. Therefore we only need to solve the 
second equation. The magnetic field intensity for region 2 will be 
completely kncwn after this is done. 

The above differential equation for i|;^^(p,z,a)) will now be 
transfonned by multiplying through by — ^ e^“^ and operating with 

yzif 

dz . 



1 






3z 



ijj (p,z,o))e^“^ 

4>2 9p2 P 9P 4>2 



= 0 



Integrating by parts twice and requiring that ip. (p,z,to) and 

9o 

■j- (p,z,o)) -»■ 0 as |z| -»■«>, we derive the differential equation 
QZ (p2 

for region 2. 






= 0 



dp 



^(p <P^^(p,a,a>)) 



= 0 



p=a,b 



( 2 . 12 ) 

(2.12.1) 
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The condition expressed in (2.12.1) follows when it is recognized that 






( 2 . 12 . 2 ) 



and use is made of the boundary condition expressed in (2.4.2) and (2.11). 

It was previously stated that the incident magnetic intensity 
does not satisfy equation (2.8.5) . In the derivation of (2.8.5) it is 
essential that the function and its derivative vanish as |z[ ->«> . 

It is obvious for slightly lossy media that does not vanish as 
z -s- -00 , therefore it cannot satisfy equation (2.8.5). 

Upon multiplication of equation (2.11) by — and operat- 

y/2^ 



ing with 



dz , the following results: 



H (p,a,co) 
Y2 



= I (hJ,) dz+^^^(p,a.0)) 



(2.12.3) 



Substituting the value of in the above integral produces: 



M(co) f 

J 



ikz . 

e dz 



(2.12.4) 



This integral does not converge. If it is assumed that the 
mediun is slightly lossy, we can write the quantity k as a complex 
number 

k = k.| + ik2 ; k.j » k2 , k^k2 > 0 

The loss term k2 will be allowed to go to zero later in the problem 
to yield results for a lossless medium. The value k is specifically 
chosen in this form to ensure convergence of the TEM magnetic intensity 
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as z ->■ +00 . This choice of k ensures convergence of(2.11.4)at the 
upper limit (z ->■ +~) . Therefore it is concluded that we should operate 

fCO 

on (2. 11) with the integral dz rather than dz . When this is 
done, the following equatiotT, which will be us'e’d later, results : 



lit (p,a,w) = + iMUO + (p,a,0)) (2,12,5) 

‘^2 v^p(a + k) '^2 



In the differential equations for the various regions of the 
problem we have made use of a quantity 




The square root is 

y =Ja - k^ (2,13) 

which is a multivalued function of the complex variable a defined on 
a two sheeted Riemann surface. The value of y must be uniquely 
defined so that only the proper branch of this double valued function 
is used. The branch cuts for this function are determined in Appendix 
B and depicted in Figure 2-3. 

The necessity of satisfying the Sommerfeld radiation condition 
requires that the asymptotic frequency domain behavior of the magnetic 
intensity for the space around the coaxial waveguide antenna (exclusive 

of regions 2 and 3 with z < 0) be that of an outward traveling or 
evanescent spherical wave. 




J 
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Figure 2-3. Branch cut for multivalued function y 
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+ z 



2 



If the value of z is such that z » p the asymptotic behavior 
specializes to 



Ha 

<Pl 

H. (p)Z,u) 
"^3 



eiMzl 



eikz 




Z 



Z -ho 



z +00 



In region 2 with z < 0 , the magnetic intensity is known [17, 
p.327] to consist of the TEM mode and higher order modes. Only 
E^^ modes are present due to the synmetry requiring no variation with 
the angle . The only field intensity components allowed are , 
Ej and . 

From equation (2.11) we have 



4). (p.z.oj) 
^2 




+ 







where , n=l,2,3,.-- 



y' is the nth ordered root of the auxiliary equation 
'on 



J (y' b) Y (y' a) = Y (y‘ b) J (y‘ a) 
o' 'on ' o'^on ' o ''on ' o'"on ' 
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C and C^„ are constants relative to the spatial variables 
p and z . J^(z) and Yy(z) are the ordinary Bessel function of the 

first and second kind, respectively of the vth order. 

In region 3 with z < 0 , the magnetic intensity is also known 
[17, p.322] to be the modes of a cylindrical waveguide. 

= j, ^on ^ 

where J\?- - = 1 /(^)^- ; n=l,2,3,.-. . is a 

constant relative to p and z and u^„ is the nth ordered root of 

on 

the auxiliary equation. 




J^(u) = 0 



Given k = k^+ ik2 , 
region 1, it is apparent 
is analytic in the strip 



and the asymptotic behavior of the field in 
from the earlier discussions that H. (p,a,w) 



-k2 < T < k2 

of the complex a plane. 

For regions 2 and 3 we must consider two cases, the first case is 

when the E modes within the region are evanescent and the second 
on 

case is when they are not. 

Case 1: ^ ^ region 2 

k < region 3 

a 
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Case 2: 


k > y' 

' on 


region 2 




u 

k >-22. 

d 


region 3 



For Case 1 the modes in regions 2 and 3 are evanescent, 
on ^ ’ 

therefore the asymptotic value of the magnetic intensity as z 

is the TEM mode portion of within region 2 and zero within region 

3. We conclude that in the event of Case 1, 



(p,a,w) 

<t>2 



and 



H (p.a.w) 
^3 



are analytic with the strip -l <2 < x < k 2 of ttie a plane. 

For Case 2, the E^^ modes are not cut off and the asymptotic 
forms behave as 



/ ^ j ^ 

-izk/l - ^ -iz(k^+ ik2)yi - \ 

e = e ^ 



z ->■ -oo 



This behavior produces a strip of analyticity that is at least given by 



-k2 < X < k2 Re(yi - (^)^) 

The symbol Re means the real part of the quantity that it precedes. 
Although at this point there is no apparent difficulty with this result, 
the inclusion of values of ka and kb such that 

ka = u„„ 

on 

and 

kb = u 

om 

where u is the mth ordered root of J (u) — J (|- u) =0 requires 
om 0 0 b ^ 

special attention. It will be shown later, that the field intensity 
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outside the waveguiding structure is zero at these frequencies. 

In light of the above, we restrict ourselves to the fundamental 
mode by imposing the conditions 

kb « 1 and ka « 1 

when we are considering a time variation other than harmonic time 
dependence. Under this restriction Case 2 does not occur, and 
the functi ons 



H 



<}>1 



(p,a,w) 






and (p,a,to) 
3 



are analytic in the strip 

-k2 < T < k^ 



of the complex a plane, see Fig. 2-4. 

When we consider harmonic time dependence, the strip of analy- 
ticity is changed to 

-k 2 < T < k 2 Re Jl - (|)^ 

- ‘^nn 

where 3 = or for region 2 or 3, respectively, n=l,2,3,---. 

The value of 3 is the magnitude of the cutoff wave vector of the 
highest order mode which is non-evanescent in region 2 or 3, as 

appropriate. Also ka u and kb ^ u 

To obtain a unique solution to the problem, it is necessary to 
determine the algebraic behavior (as a ->■ «>) of the factors in the 
completely factored and decomposed Wiener-Hopf equation [4, p,37]. 

This determination is made by considering the edge conditions [19]. 
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Air 




Figure 2-4. Strip of analyticity for (p.a.oj), 
(p,a,w) and H. (p,a,w).^ 

?>2 93 
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For the infinitesimally thin conductors utilized in this problem, the 
edge conditions are 

- 1/2 

E^(p,z,o)) ^ z ' as z->0atp = a and b 

Equation (2.5) states that 

E^(p,z,a)) « ^ 1^ [p H^(p,z,o))] 

The point of interest is that 

[p H^(p,z,o))]a: z~^!^ (2.14) 

as z -> 0^ at p = a or b . z 0^ means that z = 0 is approached 
from the direction of positive values of z . 

Development of the Wiener-Hopf Equations 

In this section the Wiener-Hopf equations will be developed. The 
solution to these equations will lead to the spatial transform of the 
desired field intensities. 

We begin by writing solutions to equations (2.9), (2.12) and 
(2.10) in that order. 



H. (p,a,u) 
<J>1 


= A{a,(ii) 


k^(yp) 


(2.15) 


\p. (p,a,o)) 
<P2 


= B(a,o)) 


I](yp) C(a, 0 )) K^(yp) 


(2.16) 


H (p,a,w) 


II 

O 


I^(yp) 


(2.17) 



The form of the solution to (2.9) expressed in (2.15) was chosen so 
that the asymptotic behavior would satisfy (2.9. 1) , thereby satisfying 
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the Sonmerfeld radiation condition. 

The solution to (2.10) expressed in (2.17) does not include the 
first order modified Bessel function of the second kind since it 
becomes infinite as p 0 a condition that violates equation (2.10.1). 
The function I,(yp) ^0 as p 0 , so H (0,a,w) = 0 and (2.10.1) 

I 9- 



is satisfied. 

The coefficients A(a,oj), B(a,u), C(a,co) and D(a,co) are un- 
known at this point. It is the task of finding these coefficients 
which leads to the Wiener-Hopf equations. Each of the coefficients 
will be expressed in terms of two heretofore unmentioned functions, 
W^(a,w) and H^(a,w) , which will be found by employing the Wiener- 
Hopf procedure. To this end we will make use of the relationships 
which exist among the magnetic intensities and their derivatives (with 
respect to p ) evaluated on the boundaries separating the three regions. 

At p = a we have 



H 



4>2 



(a,z,w) = H (a,z,w) 

(j>3 



for z > 0 



Therefore, 



(a,a,w) 

<J>2 




(a,a,w) 



(2.18) 



From the analogous condition on the magnetic intensity at p = b for 
z >1 0 , we have 



(b,a,w) = (b,a,w) 






1 



4*5 



(2.19) 



These equations relate the magnetic intensities at the bound- 



aries and will be used at a later stage in the development which leads 
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to the Wiener-Hopf equations. 

It was assumed earlier that the coaxial waveguide antenna struc- 
ture was made of a perfectly conducting material. Such material does 
not support a tangential electric field, therefore 



E^{a,z,w) = 0 
and 

E^(b,z,o)) = 0 if z < 0 

is a continuous function of p , thus E^(p,z,(o) has the 
same value at a given value of p regardless of the direction of 
approach, whether from smaller or larger values of p . 

Equation (2.5) states that 

When this is combined with the preceding discussion about the tangen- 
tial electric field, the results are 



^ [pH^^(p.z,o.)j = ^ [pH^^(p,z,«,)] 

at p = a . By a spatial Fourier transform (see equation (2.8)) of 
these equations, it is found that 



^ [pH^^(p,a,o))] 



p=a 



When we consider that E^(a,z,w) 



= 0 for z £ 0 , this equation re- 



duces to 









mf 



i 
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( 2 . 20 ) 



Note that we have used equations (2.11) and (2.12.2) to obtain (2.20). 
By very similar arguments with p=b we also find that 



These functions of a are analytic in the upper half plane given by 
T > -k2 . 

Returning to the spatial frequency domain field solutions given 
by (2.15), (2.16) and (2.17), we multiply each of the equations by p 
and take the p derivative of the resulting equations [20, p.376]. 



^ [p*! ] ♦ ^ = Yp[S(o,a.) Ij,(yp) - C(a,c«) K^(yp)] (2.22.1) 




( 2 . 21 ) 



We now define W^(a,io) and H^(a,w) as follows 



W"‘’(a,w) = -jp [pH^^(p,a,w)] 



p=a 







( 2 . 22 . 0 ) 



( 2 . 22 . 2 ) 
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By taking the last of the above equations and evaluating it at p = a , 
it is seen that 



or 




] 

3 



= W^(a,(o) 
p=a 



ya D(a,to) iQ(Ya) 



D(a,o)) 



ya IgCya) 



Note that 



[pHl ] 



dp 4> 



p=a 



since. 



E^(a,z,w) = 0 



0 



for z < 0 



(2.23) 



Recall that 

^ [pH^^(p.z,o))] = ^ [pH^^(p,z,o))] = ^ [p4'^^(p.z,a))] 

at p = a . (This may be easily seen to stem from equation 

(2.12.2) and the continuity of the tangential electric field E^). 
Since the tangential electric field is zero at p =a or b for z <. 0, 






= 0 for z < 0 



p=a,b 



and 






= 0 



P=a,b 
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Combining this result with equation (2.20) and equation (2.22.1) 
evaluated at p = a leads to the following; 



or 



W (a,o)) = Ya[B(a,o)) Iq(y 3) - C(a,o)) K^(Ya)] 



W^(a,o)) - YaB(a,o)) K(Ya) 
C(a,w) = ° 



Ya Ko(Ya) 



(2.24) 



Substitution of the values of D(a,o)) and C(a,w) given by equations 
(2.23) and (2.24) into equations (2.17) and (2.16), respectively, yields 



H (p,a,o)) ya IjYa) 



and 



W (a,o)) - B(a,o))Ya iQ(ya) 

(p,a,o)) = B(a,o)) I^(yp)-[ ^ (ya) K^(yp) 



Evaluation of the iinnedi ately preceding equations at p = a and use of 
the Wronskian relationship produces the following intermediate expres- 
sions: 

^ W'''(a,a)) I](Ya) 

4 . R(n ,,1 «*(oi,u) K,(Ya) 

= Ya yyi) ' Ya K„(Ya) 



By subtracting the lower of the two equations immediately above from 
the top equation, we find that 



ll 

JF IP 
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hT (a,a,co) - \p^ (a,a,u) + H, (a,a,o)) - \p~ (a,a,co) 



W~^(a.(o) 



B(g,(^) 

(Ya)^I„(ya) K^(Ya) Y^ K^tY^) 



The Wronskian was used to achieve some simplification in the above 
equation . 

If we evaluate (2.12.5) at p = a we find 



(a,a,o)) ■- (a,a,w) 



iM(g)) 

a(oc+k) 



But it is known from (2.18) that 



(a,a,o)) = (a,a,o)) 

q>3 $2 



Therefore, 



fit (a,a,cj) - (a,a,o)) 

<P3 $2 



iM(g)) 

a(a+k) 



and it follows that 



iM(g)) 



vTir a(a+k) *^3 



+ H" (a,a,w) - ip. (a,a,o)) = 



W^(a,w) 

v2. 



B(a,co) 



(ya) lQ(Ya)K^(Ya) Ya K^(Ya) 



Evaluating (2.22.0) at p = b we see that 



(2.25) 



H^(a,g)) 

■ Yb K^(yb) 



A(a,o)) 



(2.26) 
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Note that the boundary condition E^(b,z,w) = 0 for z 0 was 

invoked. This allows substitution of 



d 

dp 



cph: ] 



p=b 



= 0 



into equation (2.22.0). 

When (2.22.1) is evaluated at p = b it becomes 



or 



H‘*’(a,w) = Yb[B(a,0)) iQ(Yb) - C(a,w) K^(Yb)] 



C(a,w) = 



H‘^(a,w) - B(a,w) Yb iQ(Yb) 
Yb Ko(Yb) 



(2.27) 



If we replace the coefficients in (2.15) and (2.16) by the ex- 
pressions given in (2.26) and (2.27) respectively, then it is apparent 
that 

H^^(p,a,w) = - Yb^K’ifrb) *^1^^^^ 
and 

<l'^^(p*a,w) = B(a,w) I-|(yp) - My^)^^ (yp) 

By evaluation of the preceding equations at p = b and simplifying, 
we get 

^ H'’’(a,a)) K^(Yb) 

H<j^^(b,a,m) + H-^(b,a,o)) = - ' ybK~m 

+ B(a cj) H'*'(a,a)) K^(Yb) 

%(b,“,‘") ^ = Yb-TTrF) Yb l<o{Yb) ■ 
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Subtraction of the latter from the former yields 



It is seen that 

(b,a,w) - 1/;^ (b,a,w) = — iMiii) 

‘'’1 ^2 ^ b(a+k) 

by evaluating (2.12.5) at p = b and using the equality given by 
(2.19). Substitution of this result leads to 



iM(io) . u" /k \ 

^ + H. (b,a,io) 

b(a+k) ^1 



- ip~ (b,a,io) 



B(a,w) 

Yb K^(yb) 



(2.28) 



Equation (2.28) is analogous to (2.25). These two equations are the 
result of matching the tangential electric intensity (-«> < z < <») and 
the magnetic intensity (z >^0) at the common boundaries between the 
regions into which the space was earlier divided. These equations are 
also combined to form one of the Wiener-Hopf equations by eliminating 
B(a,w) between them. The resulting Wiener-Hopf equation involving 
W^(a,w) is 



iv M(gj) 
/TiF (a+k) 



[Kjj(Ya) - K^(yb)] = 




X [HT (b,a,w) - rp~ (b,a,io)] -ya K (ya)[H" (a,a,io) - i|;‘ (a,a,io)] 

^ ^3 ^2 



(2.29) 



Since C(a,w) in this problem is unique, the two va1ues--one 
given in expression (2.24) and the other given in (2.27) must be equi- 
valent. When these expressions are equated the results relate the 
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two variables W'*’(a,w) and , 



W^(a,o)) 



a K^(ya) + 

F'K (yb) -B(a,a))[j^-^ - K (ya) 

0 0 0 ' 



The terms W’(a,a)) and B(a,o)) appearing in the above expression are 
eliminated by using expressions (2.28) and (2.29), giving a Wiener- 
Hopf equation involving the function H^(a,o)) , 



iy M(g)) 
►^(a + k) 



[I^(ya) - iQ(yb)] = 



H^(a,g)) 
yb K^(yb) 



+ yb lQ(yb)[H‘^(b,a,w) - <^‘^(b,a,o))] 

- ya I^(ya)[H' (a,a,w) - i|;' (a,a,w)] (2.30) 

3 2 

In the coaxial waveguide antenna we usually have b-a « 
to ensure TEM mode excitation at microwave frequencies. For this case 
the field intensities will be rapidly varying when compared to p ,and 
the fields in the annular region will approach those of a parallel plate 
waveguide [17, p.329] of width (b-a). 

For the parallel plate waveguide it is well known from the pre- 
servation of symmetry that for a TEM mode excitation the current dis- 
tribution on the two parallel plates are equal in magnitude and oppo- 
site in direction [21, p.l26]. For excitation by E^^ mode fields, 
the currents on the two parallel plates will be equal in magnitude and 
in the same direction for odd n, and in opposite directions for even n. 




i i 
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Up to this point no approximations have been made and equa- 
tions (2.29) and (2.30) are exact. It is recognized, though, that the 
form of these equations is such that they cannot be decomposed into 
two expressions, one analytic in the upper half plane and the other 
analytic in the lower half plane with an overlapping strip in which 
both functions are analytic. This is made apparent by observing that 
the coefficients multiplying the "minus" functions on the right hand 
side of each equation are not identical. These same coefficients also 
contain a branch point and are therefore not of the form where an 
exact solution is possible [4, p.l53]. 

We note, however, that if b-a « then the field inten- 
sities approach those of a parallel plate waveguide and if TEM mode 
excitation is impressed, the currents on the cylindrical waveguide 
conductors are approximately equal in magnitude and opposite in direc- 
tion. The mathematical expression of this condition applicable to the 
structure under consideration is 

27rb[H, (b,2,o)) - H (b,z,o))] 2Tra[H (a,z,w) - H (a,z,w)] 

f -| ?2 ^3 ^2 

for z£0 . We now recall equations (2.11) and (2.7): 

H^^(p,z,w) = + t|^^^(p,z,o)) (2.11) 

^i ^ M((») ^ikz 
4) P 

Substitutions of these quantities in the above equations yields 



(2.7) 
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b[H^ (b.z.o)) - (b.z.w) - 



- a[H, (a,z,w) - ip (a,z,cj) - Mlill ] for z < 0 
93 <^2 a “ 



Under the assumption that the approximation is very good, we 
can equate the two currents and obtain 



b[H (b,z,w) - ijj. (b,z,w)] = a[H (a,z,w) - ip. (a,z,w)] 

<Pl (p2 (pj <P2 



for z < 0 



When this equation is spatially Fourier transformed the results are 
precisely what is needed to allow decomposition of the two Wiener-Hopf 
equations (2.29) and (2.30), i.e.. 




(b,a,cj) 




(b,a,o))] = a[H. (a,a,o)) 

(P3 



- \p~ (a,a,o))] 
^2 



when 



b - a « 



a+b 

2 



(2.31) 



(2.32) 



It is also noted that the expression (2.31) is analytic in the 
lower half plane given by t < l<2 or k 2 Re(yi - (3/k)^) depending 
on the excitation. It is convenient to define S"(a) such that 



S‘(a) = b[H" (b,a,w) - (b,a,w)] = a[H" (a,a,w) - (a,a,w)] 






1 



(J), 






(2.33) 



To derive decomposable Wiener-Hopf equations we substitute 



(2.33) into (2.29) and (2.30). The resulting equations are: 
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^ W'(a.a)) ^ 

v^(a+k) Y^alQ(Ya)[l<Q(Yt>) - K^Cya)] 



S”(a) 



= 0 



(2.34) 



jm + ^ s-(a) = 0 

v^(a+k) Y DK^(Yb)[I^(Yb) - I^Cys)] 



(2.35) 



It is interesting to note that if the radius of curvature is 
large we may imagine the conductors in Figure 2-1 to be cut along the 
z direction and flattened into parallel planes. Therefore it is 
reasonable to expect that equations (2.34) and (2.35) will become the 
Wiener-Hopf equations for two parallel planes in the limit as radii a 
and b approach infinity. The planes would be b-a distance apart and the 
axis of symmetry would be displaced a distance a=^b from the 
coordinate z axis. By taking the limit as a and b -> «> and using the 
asymptotic forms for the modified Bessel function, equations (2.34) 
and (2.35), respectively, become 



iM(a)) 
v^(a+ k) 



W^(a,w) 

Y^a 



Ya 



-Y- 






-] - S“(a) = 0 



sinh 



Y(b-a) 



and 



iM(co) H^(a,co) 

v^(a+ k) Y^b 



„ ^ 2 



Yb 



sinh 



1 

Y(b-a) 
~1 



S'(a) = 0 



It can also be shown from the definition of '^(a,w) and H^(a,w) 
that for a parallel plate case H^(a,w) = - W (a,w) . Therefore the 
equations (2.34) and (2.35) become identical in the limit as the radii 



approach infinity and each is the Wiener-Hopf equation for the 
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radiation from an open ended parallel plate waveguide with TEM mode 
excitation [4, p.l07 or 21, p.l28]. 

The solution of equations (2.34) and (2.35) for ^^(a,w) and 
H^(a,o)) allows a complete characterization of the magnetic intensity 
in the spatial frequency domain. The foregoing work already contains 
the equations relating A(a,aj) and D(a,o)) to W‘‘’(ct,a)) and H^(a,o)). 
These coefficients are given by (2.26) and (2.23), respectively. To 
obtain the remaining coefficients C(a,w) and D(a,w) in terms of 
'^(a,w) and H^(a,w) , we return to equation (2.22.1) and evaluate 
it at p = a and b . When these results are consi dered with (2.20), 
(2.21) and the relation 



^ (P,a,w)] 



= 0 



p=a,b 



we get 



W^(a,o)) = Ya[B(a,o))I^(Ya) - C(a,o))K^(Ya)] 



H'*’(a,a)) = Yb[B(a,w)I (Yb) - C(a,o))K (Yb)] 



The solutions to these simultaneous equations yield the coefficients 
we seek, 

H‘‘’(a,w)K^(Ya) W'’’(a,w)K^(Yb)^^ ^ 

(2.36) 



B(a,oo) = [ 



C(a,w) = [- 



yb y& 

H''’(a,w)I (ya) W‘‘’(a,co)I (yb) 



1 



yb 



ya 



^ ^Io(Yb)K^(ya) - lQ(Ya)K^(yb) ^ 



(2.37) 
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To complete the list we add: 



A(a,o)) 



-H^(a,g)) 

YbK^iyb) 






(2.26) 

(2.23) 



We return to the earlier discussion that certain values of 
ka,kb are excluded. Namely 



and 

where 

and 



ka = u 



on 



kb = u 



om 



is the nth ordered root of J (u) = 0 
on 0 



is the mth ordered root of J (u) - J (r- u) = 0 
om 0 0 b 



There is an important mathematical reason for excluding these particu- 
lar values of ka and kb that shows up in the factorization procedure 
of Appendix D. The physical reason for their exclusion is that at 
these frequencies the magnetic intensity everywhere outside the struc- 
ture is zero. 

If ka = u^^ , then the magnetic intensity in the region 
p < a , z < 0 must be zero. This value of ka corresponds to the 
cutoff wavelength for the mode of a ci rcular waveguide of radius 

a . Exactly at cutoff, the magnetic intensity for the E^^ mode in 
the ci rcular waveguide is zero [17, p.322]. 

Now considering that b-a « the structure has fields that 

behave approximately as those of a parallel plate waveguide. In this 
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limit we have shown above that the field exterior to one of the plates 

is zero. By symmetry the field that is on the exterior side of the 

other plate must also vanish. By continuity it is concluded that the 

field is zero everywhere exterior to the waveguide channel. 

For the otfier case, kb = u^ , we examine the asymptotic 

behavior of lu^l as m ^ ~ 

' om' 




with a and b of the same order of magnitude as u ->■ . Through 

the trigonometric identity 

cos X - cos y = - 2 sin(^^) sin 
it is determined that the zeros are approximately located at 

for u ^ or 
om 

The asymptotic values of u^ given above are recognized as 

the asymptotic values of the cutoff frequencies as m ~ for the 

E modes of a circular waveguide of radius ^ and the even 
om ^ 

ordered E modes of a parallel plate waveguide, with a channel 
om 

width of b-a , respectively. The second value of u^^ also cor- 
responds to the asymptotic value of the even ordered E^^ modes of a 
coaxial waveguide [22]. In the event b-a « we may approximate 
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the structure as a circular waveguide of mean radius ^ 2 ^ or a 
parallel plate waveguide, depending on the excitation. At the cutoff 
frequencies for such structures, the electromagnetic fields are iden- 
tically zero. Therefore we have excluded the use of such frequencies, 
since they invoke no response. 

Solution of the Wiener-Hopf Equations 

The two Wiener-Hopf equations (2.34) and (2.35) are solvable in 
the conventional Wiener-Hopf fashion. The various terms are factored 
into the product of a "plus" and "minus" function, where required, and 
the resulting equation is decomposed into the sum of "plus" and "minus" 
functions. [The "plus" ("minus") functions are regular in the entire 
upper (lower) half plane.] Analytic continuation arguments are then 
invoked to yield the solution. 

To preserve the continuity of the work. Appendix D was devoted 
to the factorization [23] of the kernels 

{lQ(Ya)[KQ(Yt>) - K^(Ya)]}"^ = x(a) = x'‘‘(a) x"(a) 
and 

{Kp(Yb)[Io(Yb) - iQ(Ya)]}'^ = y(a) = y^(a) y"(a) 

In this section we will use the expressions x''’(a), x"(a), y^(a), y'(a). 
Their specific functional form is found in Appendix D. 

We rewrite (2.34) and (2.35) using the newly defined symbols 
for the kernels that must be factored. 



iM(g)) ^ 
>^(a+ k) 



-45- 



[x'^(a) x"(a)] + S'(a) = 0 (2.38) 

Y a 



iM(h)) ^ H'^(g,h)) 
(a+k) 



[y (a) y (a)] + S (a) = 0 



(2.39) 



Recall that the value of 



2 

Y 



is given as 



2 

When this value of y is substituted in the above equations, we can 
reduce them to 

iM(u)) (g-k) ^ W^(g,to) x^(g) _j_ S (g)(g-k) _ q 

✓7Tr(ci+k) X (g) a(g+k) x (g) 

and 

iM((i))(g-k) ^ H~*~(g,(D) y^(g) ^ S~(g)(g-k) ^ q 

*^2Tf(g+k) y (g) b(g+k) y (g) 



Examination of the region of analyticity of the terms in the 
above equations reveals that the middle term in each equation is 
analytic in the half plane t > -kg and the last term in each equa- 
tion is analytic in the half plane t < kg . The first term in each 
equation has a pole at g = -k . If this pole were not present it 
would be analytic in the half plane t < kg . By removal of the 
residue of this pole, we obtain the completely decomposed equations: 
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W~*~(a,oo) x^(g) _ ^ kM(a)) 
a(a+k) ^ (a+k) x^(k) 

= _ ^]|( 2 . 40 ) 

X (a) v^(a+k) x"(a) x^(k) 



H^(a,g)) y^(g) _ ^ kM(g)) 

b(a+k) ^ (a+k) y^(k) 



^ _ j S~(a)(g-k) ^ iM(co) |-a-k + 2k 3 1 ( 3 . 41 ) 
y"(a) *^(a+k) y"(a) y^(k) 



In the above expression we have used the result from Appendix D, that 



and 



x^(a) = x"(-a) 
y^(a) = y"(-a) 



The terms in (2.40) and (2.41) are arranged so that the left hand side 
is analytic in the half plane T 2 > -k 2 and the right hand side is 
analytic in the half plane T 2 < k 2 or k 2 Re(Jl - (3/k)^ ) depend- 
ing on the type and frequency of the excitation. The important point 
in this discussion is that there is a strip of overlap in the regions 
of analyticity, that Is, both sides of (2.40) and (2.41) are analytic 
in a common strip 

-k 2 < T < k 2 or k 2 Re - (3/k)^ 
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We therefore conclude that each side of the expressions (2.40) 
and (2.41) represents the analytic continuation of the other side. 

Since these expressions are regular throughout the complex plane they 
are an entire or integral function. 

If we can show that each side of (2.40) and (2.41) is bounded 
as ja| -»■ 0 ° in the half plane for which that side of the expression 
is regular, then we can use Li ouvi lie's theorem to set each expression 
equal to a constant. Li ouvi lie's theorem states that the only bounded 
entire functions are constants. 

To determine the limiting behavior as [aj ~ for the terms 
in (2.40) and (2.41), we examine the edge conditions given in (2.14) 
as 



as z -»■ 0^ at p = a or b . Taking a slightly more general case, the 
Fourier transforms 






z" e^«" dz 



and 



1 

vTif 



z" e'“" dz 



determine the limiting behavior. In these integrals n < 0 to have 



absolute and uniform convergence of the improper integrals. 



It is also well known that if | f(z)| <e 



k2Z 



as z " 



then 



= lim 

(X-»-oo 



iaF^(a) 



(lim ia)( lim F^(a)) 
cx-^- «> a-^ 



lim f(z) 
z-vO+ 
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kpZ 

for T > -l <2 where a = a+ix . Similarly, if |f(z)| < e ^ as 
z -V -00 then 

lim f(z) = (lim ia) (lim F"(a)) 

Z->0*" (X-)--oo 0^->_oo 



for T < k 2 . Thus the limiting behavior as (a( ->-00 in the region 
of regularity of the spatial transforms is determined by the behavior 
of the space domain function in the limit as the origin is approached. 
Using the ganma function [18, p.l83] 



r(x) 





6 



Re X > 0 



and contour integration, we find 



and 



00 




eiaz 



dz 



+ 



1 

v^ZtT 



r(n+i) 

n+1 

a 



i f(n+l) 
e 



1 



0 



laz 



dz 



1 r(n+l) '' 



with -1 < n < 0 in both integrals. The lower extremity on the 
range of n is set at -1 to ensure the existence of the gamma 
function. 

From the edge condition we have n=- ^ and determine 

..+ / ^ -1/2 
W (a,a); a 

and 

H‘^(a,w) as ja] with x > -k 2 . 
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From the definition of S~(a), it can be seen that it is the 



spatial transform on negative values of z of the current normal to 
the edge at z = 0 on either the inner or outer conductor, exclusive 
of that current contributed by the incident field. The magnetic in- 
tensity is continuous at all values of p for z = 0 , thus the 
current in each conductor must approach zero as z -»■ 0~ • This is 
apparent from the fact that the current in each conductor is propor- 
tional to the discontinuity in the magnetic intensity across the con- 
ductor. The discontinuity disappears at z = 0 . Because the total 
current approaches zero as z ^ 0~ , the quantity which transforms to 
S~(a) must approach that of a TEM mode field intensity traveling in 
the minus z direction. This quantity should just cancel the incident 
mode field intensity at z = 0 . Thus 




=> M((jo) e~^^^ as z -»■ 0 



The transform of this quantity is 



0 




From this result we conclude that S~(a) ^ ^s ja| -> «> with 

T < +k2 . 
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Collecting some results from Appendix D, we write 




1/2 



and 




Also, 



x"(a) 



and 




An examination of (2.40) and (2.41) in the light of the 
asymptotic behavior of the terms shows that each side of these ex- 
pressions goes to zero as |a| in their half plane of regularity. 
Liouville's theorem is therefore applicable and each side of (2.40) 
and (2.41) must be identically zero. The quantities of interest are 




and 




Therefore 




(2.44) 




(2.45) 
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and the coefficients A(a,w), B(a,to), C(a,w) and D(a,w) have been 
determined, see equations (2.23), (2.26), (2.36), and (2.37). The 



spatial domain field intensities are now written: 

k M(g)) 

YKQ(yb) y'^(k) y'^(a) 



H^^(p,a,o)) = - 



(2.46) 



H^3(p,a,(o) = ijl 



- ^ /2 k M(g)) M 



I-|(yp) 
+7ri +"/ 



YlQ(ya) X (k) X (a) 



(2.47) 



and 



i|»^^(p,a,w) = *-yYb)yYa)^- iQ(ya)yYb)^ 



^I](yp)KQ(ya) + K^(yp)I^(ya) I^(yp)K^(yb) +K^(yp)I^(yb) ^ 

X I - . ~ - j } 



y^(k) y"^(a) 



x^(k) x^a) 



(2.48) 
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3. Frequency Domain Description of Fields Interior to the Coaxial 
Waveguide 

In this chapter we shall derive the frequency domain fields 
interior to the coaxial waveguide. We will completely characterize 
the TEM mode as well as the (q=l ,2,3, • • • ) mode fields. As pre- 

viously stated no other modes are allowed due to the symmetry of the 
structure . 

Let us begin by recalling equation (2.11) 



H. (p,z,w) 

^2 



+ (p,z, 0 )) 

(f>2 



The incident field is known to be 

^i _ M(g)) e^*^^ 

<j) p 

Our task is to find ip (p,z,w) . From the definition of the spatial 
transform we write 



ip, (p,z,w) = — ^ ip, (p,a,u) e da (3.1) 

<J>2 J ^2 

.oo 



Before attempting this integration we will examine the integrand 
and all singularities will be located. To do this we make use of the 
relations [20, p.375] between the modified and ordinary Bessel func- 
tions that follow: 



1 



1 



- y vin ^ Tri 

I^(z) = e J^(z e'^ ) ; 



1 

K^(z) = ^ Tii e^ 



1 n 

VTri / 1 \ T 1 

Hy^ze^ ); 



(-IT < arg z £ |-) 
(-11 < arg z £ |-) 
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We also make use of the analytic continuations [24, p.80] 



I (z e"^M = I (z) 

V V 



K (z e"’) = K (z) - i ^!L«£ i (z) 

v' ' v' ' sin vir v' ' 



By these formulas it can be shown that 



Io('Vb)Ko(ya) - lQ(Ya)K^(Yb) = - |£jQ(Y'b)YQ(Y' a) - J^(y' a)Y^(Y'b)] 



where 




It is clear that expression (2.48) does not have a branch point 
at a = +k , since we may replace y by -y and not change the value 
of the expression. However, (2.48) has a branch singularity at 
a = -k due to the presence of x^(a) and y^(a) . From the above, 
we conclude that the only singularities of ip. (p,a, 03 ) which exist in 
the upper half plane are poles, located at a = k and a = 

(q=l ,2,3, • • • ) with 



ot. 









where y‘ 

' oq 



is the qth ordered root of 



JoW'b)Yo(Y'a) - Oo(Y'a)YQ(y'b) = 0 (3.2) 

The determined by (3.2) are the eigenvalues determining the cut 

off frequencies of the modes of a coaxial waveguide [22], 

To evaluate the integral (3.1), we integrate along tiie closed 
contour in the upper half plane as shown in Figure 3-1. On this 
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Figure 3-1. Integration Contour 
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contour the integral converges for z < 0 and the contribution along 
the semicircular path at |a| ~ is easily seen to be zero by re- 

placing the modified Bessel functions in ip, (p,a,w) by their 
asymptotic values and using Jordan's leima to evaluate the resultant 
integral over the portion of the contour at ja| -> <». 

By the residue theorem, we have 

Ip, (p,z,w) = 2-iri I residues of ip. (p,a,w)e“^“^} 

at the poles interior to contour 
of Figure 3-1 (3.3) 



The residue of the pole at k is 



i,\rikM(co)e nr 1 i 

im(a-k)[ ^ (Yb)K (ya) - I (ya)K (yb)^ 



a->- k 



O'' ' 0 



0 '' ' 0 



I-|(yp)Ko(Ya) +l<i(yp)Io(ya) (yp)K^(yb) + (yp)I^(yb) 

^ + 77 ^ 



y (k) y (a) 



x^(k) x^(a) 



} 



_ ik M(g))e~''‘'^^ iim[»CT] 
"" a->k^d¥ 



Cy^(k)]'2(^ iln ya+^) - [x-^(k)]-2{- iln yb + 



JC£ 
)ln b/a 






-ikz 



(3.4) 



where limiting forms of the Bessel functions have been used to charac- 
terize their behavior for small arguments. 
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The residue of the pole at a Is 

q 



lim (a-a ) 



ikM(g)) e 



-laz 



a-»-a Tnr{- 



^I^(Yp)K^(Ya) + K-j{Yp)JQ(Ya) I-](YP)l<Q(Yb) + K^(yp) JQ(Yb) ^ 
I , + , . + '. . + . . } 



y^(k) y'^<a) 



x^(k) x^(a) 



Using the L'H6pital rule this becomes 
-iz 

lim [. Mw)e 2 

a^aq ^ 



X< 



^ ^Ji(y'p) YQ(Y'a) - Y^(Y'p)JQ(Y'a) Ji(Y'p)YQ(Y't>) - Y^(Y'p)JQ(Y't>) 
y^(k) y^(g) x~^(k) x~*~(a) \ 



|^[Y'(-f)(J,(Y'b)Y^(Y'a) - J^(Y-a)Y^(Y'b))] 



(3.5) 



Observe that y' = /?r7 and y' = v/k^- ct^ . By invoking relation 

oq q 

(3.2) the derivative in the denominator of the above relation is readily 
reduced to 

^ (y'[- f][0„(Y'b)V„(Y'a) - J„(Y'a)V„(Y'b)]} 

= - f Y'[|^J„(Y'b)Y„(Y'a) - 0^(Y'a)V^(Y'b))] 



= - jY'[^Kb(0,(Y'b)V„(Y'a) - J„(Y'a)v„(Y'b)) 
+ a(jQ(Y'b)Y-|(Y'a) - J-| (y' a)YQ(Y'b) )] 
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Putting this value in (3.5) yields the residue 



-i z (y‘ 

1 k M(g)) e ^ 

(y‘ 









y*(k)y*i/k^- (Ypq)^) x'^(k)x'''(yi<^- 



(3.6) 

J (y* b) Y (y* 5) 

The relation used to simplify (3.6). This 

o^'oq ^ o^'oq ' 

relationship follows from (3.2). 

(p.z.cj) for z£0 is found by substituting (3.4) and (3.6) 
into (3.3). The result is 



oo -1 z/k^-J^' 

ijj^^(p,z,u)) = -M(co) M(w)R^f^(p)e 

z < 0 



■ ST¥7? 



(3.7) 
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^ ^ 



1 J (y' b) 

y^(k)y~*~(s/k^- (y^)^ ) x~*~(k) x~*~( v/k^-(Y^^)^) I 

b(J,(Y' b)Y(y' a)-J (y' a)Y,(y' b))+a(J (y' b)Y,(y' a)-J,(y' a) Y (V b))( 
r^oq ' 0 *oq o'^oq ' l''oq ' ' o 'oq * P''oq ' r^oq ' o'^oq 

(3.8) 

f^(p) = Jl(v;qp)Y„(Y;^a) - Y,(v;qP)Jo(Y;K;a) (3-5) 



The magnetic intensity for the interior of the feedline is seen to be 



H. (p.z.cj) = 6 — + if) (p,z,w) , Z£0 (3.10) 

Y2 ^ ^2 
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4. The Radiated Electromagnetic Field 

To find the magnetic intensity in region 1, we must take the 
inverse spatial Fourier transform of the spatial frequency domain mag- 
netic intensity found in Chapter 2. The integral to be evaluated is 



for p^b . This equation was obtained from equation (2.46). 

Examination of the integrand of (4.1) using the analytic con- 
tinuations of the Bessel functions and the value of y on both sides 
of the branch cut reveals that the integrand has branch points at 
a=±k. Any attempt to evaluate the above integral in closed form by 
contour integration would have to account for integration along a path 
alongside the branch cut. If the contour is closed in the upper half 
plane (for z £ 0) then the contribution from the path alongside the 
branch cut in the upper half plane must be calculated. If the contour 
is closed in the lower half plane (z >. 0) , then the contribution from 
the path alongside the branch cut in the lower half plane must be 
calculated. In the special case of z = 0 the integration around the 
branch cut in the lower half plane resulting from closing the contour 
in a negative sense equals the integration around the branch cut in 
the upper half plane which results from closing the contour in a posi- 
tive sense. There are no poles located inside the contours under 
consideration . 




y k M(cj) K^(Yp)e"^“^ 



(4.1) 



_i YKQ(yb) y"^(k) y'^(a) 
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The integrations around the branch cuts are quite complicated 
and are not evaluated in closed form. Since our primary interest is 
in the radiated fields (i.e., those with 1/r behavior), it is sufficient 
for our purposes to approximate the value of the integral (4.1) for 

-► OO 5 

r = 

In the limit as r ~ , we replace the Bessel function 
dependent on p by its asymptotic value for large arguments [20, 
p.378]. Using Figure 4-1, we readily have 



p = r cos X 



and 



z = r sin X 



where |x| <. ir/2 . Using the asymptotic value 



k^(yp) 

the integral (4.1) becomes 



♦i ^ y+(k) i 



'\j 



TT „-YP 



2yP 



g-(Yr COSX+ ictr sin x) 
yKq(yU )y"^(a) [ /ytToTx ] 



da 



as r ^ ~ , p i b 



(4.2) 



The variable 6 has been introduced at this point since this 
is the variable in which we wish to have the final expressions. It 
is obvious from Figure 4-1 that 



IT 



7 ■ 



X 



e 
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Figure 4-1. Coordinate Axes 
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Define a change of variable by introducing the transformation 

a = k sin z (4.3) 

z is the complex variable x+iy . Under this change, the function Y 
becomes 

Y = \/k^sin^z - = -ik cos z (4.3.1) 

It is necessary to specify precisely the value of z , since it 
is given by the multivalued function 

z = arc sin(^) 

which has infinitely many branches. The precise specification will 
allow a one to one mapping from the complex a plane to the complex 
z plane. Since 



we may write the equation 

e^^ ^ - 2i ^ e^ ^ - 1 = 0 

This quadratic is readily solved for z 

z = -i zn[1 f + [1- 

In this solution, the principal branch of the logarithm function has 
been used and the plus sign has been used with the square root function 
[25]. By the branch previously defined for y » the above equation 
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simplifies to 



2 = -i tn[i 



(4.4) 



The equation just developed is useful in determining how the two- 
sheeted a plane maps into the z plane. 

Under the transformation (4.4) the numbered regions of Figure 
4-2 map into the correspondingly nunbered regions of Figure 4-3. The 
path of integration for the integral of (4.2), depicted in Figure 4-2, 
maps into the path shown in Figure 4-3. 

As examples consider two points on the path of integration in 
the a plane 

1) a wi th I t| 0^ 

2) with |t|->0 



For case 1: 

Y = lim [(-|a| + i | t| )^- 

a 

M->o+ 



= lim [|a|(l 

Q -^OO 

|t|->0+ 



2|a 



t) - ih|] 



and 



2 = -i lim iln[i ^~l^ l _-j -i-jm [^,n 

|a|^ \o\^ 

|t|-K)+ 



-ik 

2R 



-] = 



-2. + ioo 



For case 2: 



2 

Y = lim [(|a| -i|T|)^- = lim [|a|(l - — — ^) -i|T|] 

|a|-»<“ 2|a| 

|t|->0" |t|^ 0" 



and 
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^3 Inri y >0 

I I Im y < 0 




Figure 4-2. Integration path in a plane 
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Figure 4-3. Integration path in z plane 
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z 



-1 lim 

lal-^ 

IxhO- 



-IX 






-i lim iln[i = j 

laH 2 






Upon substituting (4.3) and (4.3.1) into (4.2), the integral 
under examination becomes 



(r,e,co) 



M(a3) k 
v^y^(k) 



oo 

f exp(ikr(cos y cos z - sin x sin z)) 

■* y^(k sin z)K^(-ikb cos z) /-irk cos x cos z 



-oo 



(4.5) 



The integral is in the form amenable to the saddle point method of 
integration, since r . 

Define the z dependent quantity in the exponent as 

g(z) = i(cos X cos z - sin x sin z) 

= i cos(z + x) 

It follows that 

^ g(z) = -i sin(z + x) 

The saddle points z^ are the points at which g(z) is a maximum. 
Therefore from the above we know that 



Zg = htt - X 



n=0,±2,±4 ••• 
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Con centrating on the point z = -x » we note that g{z ) = i . 

s s 

Since a steepest descent path is a path along which the imaginary part 
of g(z) is constant, the steepest descent path through = -x is 
given by 

Im g(z) = 1 

Since g(z) = i cos(x+iy + x) = i(cos(x+x)ch y - i sin(x+x)sh y) , 
the steepest descent path through z^ = -x is given by 

cos(x+x)ch y = 1 

This path is shown in Figure 4-4. We are assured that this is a path 
of steepest descent rather than one of steepest ascent, since the 
Re g(z) monotoni cal ly decreases along this path as we move away from 
the saddle point. 

By the residue theorem, we defonn the path of integration of 
Figure 4-3 into the steepest descent path. The major contributions to 
an integral along a path of steepest descent come from points in the 
neighborhood of the saddle point. The slowly varying terms of the 
integral are replaced by their value at the saddle point and the func- 
tion g(z) is replaced by the leading terms of its Taylor series about 
the saddle point. 

g(z) - i - 2 ^z+x)^ 



M(0))k dz 

(r,9,o)) == ^ ^ 

^ y (k)y (-k sin x)/-i cosyK (-ikb cos x) 



(4.5) 




■ 




1 



Pi 



I 
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Figure 4-4. Steepest descent path 
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Along the steepest descent path 



^ ^ ~ cos (x+x) 



and 




+ , t anj x ,tx) 

cos(x+x)s 



nu+y; 

(x+x)sh y 



As we approach the saddle point in the manner x->-x and 
y 0^ , the slope is 




As we approach the saddle point in the manner x -+-x and 
y 0" , the slope is 



Therefore the angle that the SOP makes with the real axis at the 
saddlepoint is either 

T - 4 

Examination of Figure 4-4 shows that at the saddle point the SOP and 
the negative real axis form the angle 3 tt/ 4 in the second quadrant 
and -tt/ 4 in the third quadrant. This observation permits us to 




express 




z + X = u 



I u| e ' ; u e SOP and y > 0 

|u| e"^ ’ u e SOP and y 0 
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It is also observed that 



e~^ at the saddle point 



Since equal contributions to the integral of (4.5) are received as we 
move away from the saddle point in the two directions of increasing u, 
we write the integral which appears in (4.5) as 



-i I 



” _ kri |2 

exp(-i ^(z+x)^)dz ^ Z I e e ^d|u| 

SbP 0 



(4.6) 



Using the well known integral 



.2 2 r- 

e'3 a dx = ^ 
e ax 2a 



the value of (4.6) is seen to be 






2 (e 



^lu'l 



- i ^ 



') d|ul = e’' 



Vkr 



(4.7) 



Equation (4.5) becomes 



H, (r,6,w) = M(a)) 



j'kr 






1 



l-r-T 



y (k)y (-k sin x)l<Q(-ikb cos x)cos x 



} (4.8) 



To simplify the above expression we use the identity relating modified 
Bessel functions with ordinary Bessel functions: 



K^(-ikb cos x) = i |-H^^^(kb cos x) 








I 








From Appendix D, the expression 
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y(a) =[i f H<’>(Kb)[J^(Kb) - Jjj(Ka)]T 



is formed, where K = +iy . In expression (4.8) we have 



i cos x) = i sin 6) 

If in the expression y(a) , K is set equal to k sin 0 and a is 

2 2 

set equal to k cos 6 , which follows from the fact that y =“K = 

2 2 

a - k^ , then 

y(k cos e) = y^(k cos e) y~(k cos e) 

= ji 2 'HQ^^(kb sin 0)[JQ(kb sin e)-J^(ka sin e)]J 

When this equation is solved for i ^H^^^(kb sin 9) and the result 
substituted into (4.8), the far zone magnetic intensity is seen to be 



H . ( r , 0 ybi) 
<}>1 



= M(w) 



ikr ^ ykb sin 6) - J^(ka sin cos e) , 

y^k) 

(4.10) 



sin e 



The relation (4.10) has the restriction that p ^b . In the 

region z < 0 the equation (4.10) is only valid in the region p >ib; 

however, for z > 0 , the H. can be analytically continued from 

<?>1 

region 1 to 2 and from region 2 to 3. Therefore the restriction re- 
quiring p >. b can be removed in the region z >_ 0 and equation 
(4.10) becomes an expression for the magnetic intensity for all values 
of p in the far zone (r ->-«>) exterior to the feedline. 
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It is recognized that equation (4.10), except for the factor 



y^(k cos 9) 
y*(k) 



, exactly equals the results obtained from assuming that 



the incident field is the total aperture field and that no contribution 
is made to the radiated fields by currents which flow on the waveguide 
structure . 



waveguide antenna first obtained in this work modifies the field pat- 
tern obtained by the conventional solution to the problem of radiation 
from this type of antenna [19,p.302]. It will be seen later that this 
factor has a significant effect at high frequencies and modifies the 
field pattern to direct the radiation in a more forward direction. 

By equation (2.11) we expressed the magnetic intensity in 
region 2 for all z as the sum of the incident field and another term 
which represented the reflected and higher order mode field intensities. 
By the radiation condition, we know that the incident field cannot 
exist in region 2 for z > 0 . 

To show that the incident field intensity is canceled in region 
2 for z > 0 , we evaluate the integral 



by using the contour of Figure 4-5, ip. (p,a,w) is given by (2.48). 

4>2 

The value of the contour integral is calculated by the residue theorem. 



The new field factor contribution 




to the coaxi al 



oo 




(4.11) 



The value of the integral (4.11) is the negative of the contribution from 
the path around the branch cut and -2iri times the summation of the 
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Figure 4-5. Integration contour 
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residues at the poles of the integrand located on the interior of the 
contour c. The contribution along the semicircular arc at ja| ^ is 
zero. 

The integral (4.11) is not evaluated in closed form due to the 
complex integration on the path alongside the branch cut. However, we 
shall show that the contribution to the value of the integral from the 
circular path around the branch point (k) is just enough to cancel the 
incident field intensity. 

"i fi 

On the circular path around the branch point a = -k + ce , 
e 0 and 0 £ 0 £ 2ti . The integral of interest becomes 



lim 
e-^- 0 



J '^2 




-k 




lim 

e->0 





We have made use of the limiting forms of Bessel functions for small 
arguments and the following: 
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x^(-k) = x’(k) ; x^(k) x”(k) = x(k) 

y'^(-k) = y'(k) ; y'^(k) y"(k) = y(k) 

To see that this is just enough to cancel the incident magnetic inten- 

sity, note that 

1 1 _ b 

ylFT' ^ ? 



which is easily found from the definition of y(k) and x(k) given 
in Chapter 2 and the limiting forms of Bessel functions for small 
argunents. 





i 
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5. Time Behavior of the Radiated and Reflected Electromagnetic Fields 

The time behavior of the field intensities is found by comput- 
ing the inverse time transform of the frequency domain expressions. To 
find the time behavior of the radiated magnetic intensity emanating 
from the waveguide antenna, the inverse transform of equation (4.10) 
must be obtained. Because of the complicated functional behavior of 

the factor y cos e) computation of the integral 

y (k) 

H (r,e,t) = 

X e-'“t do (5.1) 

y (k) 

is extremely difficult. 

For harmonic time dependence the exciting voltage V(t) is 
iw t 

Re(A e ° ) (see equation (2,7)). Then 



M((o) 



ikr J (kb sin e)-J (ka sin 0) 



[- 



sin 0 



-] 



H(») = «(»-»„)] 

where 6(x) is the Dirac delta function. Therefore, the case of 
harmonic time dependence (5.1) reduces to 



H (r,0,t) = Re 






y"(^cose) -'K(t-r/c)] 



-^^0 lu b/a y^{^) 



r sin 0 



(X) b 0 ) a 1 

[J^(-|-sin 0) - Jg(-|-sin 0)]J 



X 



(5.2) 







- 1 
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The value of k is w/c as given by (2.4.3). 

+ “o 

The factors of equation (5.2) except ^ are well 

/<^) 

understood, and are not in need of further simplification. 

From Appendix D we transpose the function 

k_ 



+ , , » oo 

y (k cos e) ^ 
y'^(k) m=l 



r 1 +r 



m 



1 + 



k cos 9 

'Kl ■' 



, (5|^l-cos e) i^sin e 



exp 




K^^^u) in 



1 + 



1 + 



k cos 6 



- - u^- 




,a+b 



k(fg)(cos e- l)(l-C+i 



exp^k(j^^cos e[-)ln(k(|^)cos e) + ^^«,n(2 cos e) 



(5.3) 



In this formula we have used the fact that if a = k cos 6 , then 



«_/Y+(a-k)\ _ „„r-ik sin 6 + k(cos 0 - l)i _ «..re~^®-11_ ,• 
Ta ' -kT * - sin e - k(cos e - 1)J - “"iTT^e J ‘ 



0 



using the principal branch of the logarithm function. Additionally, 






3 S >] 

/(ub)+Y:(ub) 



ro|=4 
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It is also known from Appendix D that y^(a) ^ as 

a ^ . From this asymptotic behavior we see for very high frequen- 

cies (X « b-a) that 

lim |cos (5.4) 

k * » y (k) 



In Chapter 2, there was imposed a restriction that kb,ka « 1 

for non-harmonic time dependence. In keeping with this restriction, 

equation (5.4) should only be used for the case of harmonic time 

dependent excitation. For values of ka,kb « 1 the value of 
y^(k cos 0) 

— is more complicated but instructive to determine. To 

determine the value of this field factor in the low frequency limit 
(X » b-a) we use the series expansion of the infinite product 



ii|^i-cos e) 



1 1 + 


k 


\ ^ i 


1 a 1 
' m' 


1 + k 


cos 6 


y i\ 


a 1 
m' 



= 1 - ik(^|^)(l - cos 9) y [ \.L ^ 



' ] * 0(ka)^ 



The infinite series is the difference of two divergent series but is 
itself convergent [26, p.33]. This is apparent since 
1 J_ ^ 1 _ J_ = _L 

miT 



a. 



a+b 



m' 



(m - ^)n 



J_ 

mir 



4ivm 



1 

and the infinite series ^ 

m=l m'^ 



as m “ 



is known to be absolutely 
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con vergent. The following integral 



+ K^^^(u) £n 



6M) 



1 + 



vgT7 

1 ^ k cos 6 



du 



also appears in ^ . This integral will now be shown to be 

y (k) 

convergent and its approximate value will be determined in the limit 



kb,ka « 1 . 

Upon substitution of the value of K^^^(u) and simplifying, 
the above integral becomes 



«(2), 






6^ 



j2(ub)- Y^(ub) 



■)] 



fi + J _.-. cos e _ 1 

cos 6 +/l - 



du 



Introducing a change of variable, ub = z , we get 



1 

TT 



6 ->0 



1 - 



7 TZ ^ 



j2(z)tv^(z)^. 



£n 



1 + 



1 - cos 0 



cos 0 + 






dz 

(5.5) 



To examine the behavior of (5.5) for small values of z , we 



consi der 



TT j . 7 TZ 

cS>0 



1 



7 TZ T . 4 , 2 

1 + ^ £n z 

TT 



] £n 



1 + 



1 - cos 0 



cos 






dz 
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This integral converges as 
Now 



e 

r 






1 



7T 



)] dz 



converges . 



e 




0 




dz 

z J,n^(z) 



1 



7T 



4 an z 



e 

0 



This quantity is seen to be convergent for all e ^ 1 . We can also 
remove the principal value sign and change the lower limit to zero 
without affecting the value of the integration. 

At the other extreme for large values of z , we consider 









1 + 



1 - COS 0 



dz 



0 



cos 0- - 1 



for M » 1 . The asymptotic value 

i(z - 1 ) 






^ [1 + o4)] for z -V ~ 



is now substi tuted wi th the result 






-Lp-] )ln[ 

+ o(-!-n L 



1 + 



1 - cos 0 



[l+0(^)] 



1 

TT 



M 



0(|) an[l + 



cos 0 
1 - cos 0 






.] 



dz 



cos 



e 






)'- T 



dz 
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The logarithmic function iln(1+x) can be expanded 

2 3 

£n( 1+x) = x- -^+^-*-- , |x|<.l; xyf-l 

When this expansion is considered along with the fact that 



1 + 



1 - cos 6 



= 5 > 1 + 



cos 0 - i - 1 



as 



Z 



«> (with z real and |A| finite), then 



i f 0(1) ln[ ^ ^ Q ] dz - 

IT 1 'z j 

M cos e - i^) - 1 



oo 

1 f 



IT 



M 



0(1) (|) 



dz 



This integral converges. 

The integrand of (5.5) has an integrable singularity at z = kb 
when 0 = it/ 2 . By the above we conclude that, for all values of 0 
and k the integral of (5.5) is convergent. 

y^(k cos 0) 

Since our interest is in the value of — y~(lD 

ka,kb « 1 , we next endeavor to find the value of (5.5) in the limit 
as kb ->■ 0 . We will use this limiting value as an approximation for 
the integral appearing in equation (5.3). 
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lim - 
kb ^0 ^ 



[1 



4 <,z 



o;(z)*y‘(z) 



)] 111(1 + 



1 - COS 0 



cos e + '^1 - 



dz 



kb 

lim — [ [1 - — (-9 

kb^o""! J^(z) +Y^(z) 

0 '0 



0 



- — 2 )] + 

+ Y ^( z ) V 



1 - COS 0 



COS0 + yi - (|^)^ 



dz 



+ lim i 



[1-§t( o y ^ )]^n(l + 1 ~ Q 



kb.O"^, J2(z)n^(z) 



cos 0 - i ,/(|g-)^- 1 



dz 



Note that the quantity 



1 - cos 0 



1 i m [8,n ( 1 + 

kb ^ 0 - . //TT2 7 

cos 0 -1 /(j^) - 1 



)] 



From the above work we conclude that as kb goes to zero, the integral 
(5.5) goes to zero. Thus, for low frequencies, we approximate the 
function 



exp 



00 

f 



6-i'^ 



K^^^(u) )ln 



1 + 



y ?:7 



1 + 



k cos 0 



du 



appearing in equation (5.3) by unity. 

If in (5.3) we make use of the limit that x £n x = x 

then the field factor y 0) approximated as 

y (k) 



as 



0 , 
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(l . - cos ej J 0(ka)^) 

y (k) m=l |a„| -TT- 



'm> 2 



exp ^ i sin 0 + ik(|^)[C(l - cos 0)]- ^^cos 0 £n(2 cos 0) - £n 2]j 



for ka,kb « 1 



(5.6) 



The time behavior of the far zone radiated magnetic intensity 

for a non-harmonic time dependent incident field intensity with 

ka,kb « 1 is given by (5.1) when equation (5.6) is used for 

y (k cos 0 _ ) _ ^ 2one the value of r in equation (5.1) is 

^ a+b 

very much larger than b or -y . To the approximation used, the 
exponential function in (5.6) is replaced by one and for that reason 
it does not contribute to the phase factors in integral (5.1). In 
addition the term of (5.6) containing the infinite sum can be dropped. 

This is apparent since ka,kb « 1 and the Fourier transform of this term 
will be of order 1/c , which will make its contribution very much 
smaller than the contribution from the remaining term. In view of this 



and 



y^(k cos 0 ) 






M(g))e 



1 



-i (wt- kr) 



r sin 0 



[J^(kb sin 0)- J^(ka sin 0)] dw 



for ka,kb « 1 



(5.7) 
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When the incident magnetic intensity is harmonically time 
dependent, then the value of (5.7) is easily obtained. In this case 



M(cj) = Re*/^ 



— — - — 6(o)-(i) ) 

^^0 an b/a ° 



and for low frequencies (A » b-a) the radiated fields are 



fe~ cos(oj (t- -)) w b w a 



EQ(r,0,t) = 






(5.8) 



For non-harmonic time dependence the radiated field can be 
found by direct integration of (5.7) or by a time domain convolution 
The integrand of equation (5.7) can be considered as the product of 
the Fourier transforms of two time domain functions. Therefore this 
integral can be evaluated by the convolution theorem and the impulse 
response is readily identified from the result. 

In applying the convolution theorem, use is made of the fol- 



; 0 £ 1 1| < z 

; |t| = z 

; otherwise (5,8.1) 

Upon evaluation of equation (5.7) by the convolution theorem we get 



lowing well known integral [27]: 



^ /2/it 

x/?T7 



[ J (az)e~^^^ da ={ 
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the result which follows; 



f 



“iFTnre I * 

— OO 

where f(x) is 



(5.9) 



sin e)^- 



f(T) = 



v/(^ sin e)^- y& sin 6)^- 



; T < - — sin 9 
c 



; - ^ sin 6 < T < - ^ sin 0 



;; - sin 6 £ T _< ^ sin 6 



sin 0)^- 



; sin 6 £ T £ ^ sin 6 

L. C 



; T > — sin 6 
c 



(5.9.1) 



and 



f = t - 



The response to an impulse of voltage excitation is therefore 



Fo 

H^(r,0,t) J Trr sin 0 In b/a 



(5.10) 



From the character of the impulse response we see that at an 



observation point in the far zone, the radiation lasts a total time of 
(2b/c)sin 0 . By ignoring for the moment the frequency restriction of 



1 



c 




I 








hi 
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ka,kb « 1 and allowing an excitation with zero rise time, a step 
function of voltage , then by integration of the impulse response we 
obtain the step response. 



H^(r,0,t) 




Trr sin 0 Jin — 
d 



X 



IT X T C 

7 b Sin 0 



T C TC 

b Sin 0 ~ a sin 0 



TT X ^,-„-l T c 
- y + sin v-v.y --- Q 
L b sin 0 






where x = t - — . 

c 



; T < - ^ sin 0 
; - ^ sin 0 £ T £ - sin 0 

; - sin 0 £ T £ sin 0 

; sin 0 £ T £ ^ sin 0 

; T > ^ sin 0 
— c 

(5.11) 



It is advantageous to have the step response for numerical cal- 
culations since its magnitude is finite at all points and the far zone 
fields can be found by the convolution of the time derivative of the 
excitation and the step response. This procedure does not violate the 
low frequency restriction. 

Because of the low frequency restriction the input excitation 
should not contain any high frequency components. To render the high 
frequency components in the excitation negligible, the rise time of 
the excitation pulses must be long compared with the time b/c . If 








m 
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the input is of the form [1 - u(t) then the time constant T 

must be chosen so that T » b/c . This will preclude any violation of 
the low frequency restriction. 

For such inputs, equation (5.9) contains convolution integrals 
of the type 

J - J -- ~ --- dx, with - — sin e £ T £ - sin e 

- ^ sin e '/(^sin e)^- 

and T < t' 



- It-Tl 

(1 - e ^ ) dx ^ 



.M 

( 1 - e dx 



.b , .b e)2- 

c c 



|t'|/T 



t' 



X dx 



--sine A- sin e)^- 






dr 



sin 

c 



e ./(-sin e)^-T^ 



e-|t'i/T/-r 5 = 

■^--p-y(|sin e)2-(t')' 



0 < |t'| £ I sin e 



Since T » b/c we neglect the second term and conclude that 
the response to an input of finite rise time with time constant 
T » b/c approximately equals the step response multiplied by 

(1 - e'l^l ! . 
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Exci tati on 




Figure 5-1. Radi ated magneti c intensity of the coaxial waveguide 
antenna for (a) step function and (b) exponential excitation, 

/s K 

Ordinate: H. ^ /i. r £n - sin 0) (volts) 

(f) 9 0 0 d 

Time constant for (b) is cT/(b sin 0) = .1 
Antenna Dimensions: a .8b 
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A basic assumption used in Chapter 2 is that the total current 
in the inner conductor is equal and opposite to that in the outer con- 
ductor. This assumption implies that the contribution to the radiated 
fields from all parts of the structure outside of the aperture are neg- 
ligible. Examination of the impulse response shows that the radiation 
appears to originate from the edges of the open ended structure 
located in the aperture plane. The outer edge gives a positive con- 
tribution and the inner edge gives a negative contribution. The 
radiation starts when the incident current reaches the edge and lasts 
a time proportional to the time it takes light to traverse the diameter 
of the conductor. The difference in sign of the contributions is due 
to the difference in sign of the current on the conductors. 

We now turn to the problem of determining the time behavior of 
the field intensities for a_<p_<b; z^O . The magnetic intensity 
in this region is given by 



H. (p.z.t) = 



/B 



M(m) 



etkz „^-1kz 



-izA^- (y')^ 

- I R f (p) e 
q=l ^ ^ 

-iw t 



dco (5.12) 



For harmonic time dependent excitation voltage Re(Ae ° ) , 



H(u,) = Re[/S 

0 in 



u, ' O' 



and (5.12) becomes: 



I 



^ I 



o 



i 

I 







-90- 




R(o)o) 



(5.13) 



where R(a)^), Rq(<*JQ) fq(p) the functions (3.7), (3.8) and 

(3.9) evaluated at o) = w , if applicable. oj should not corres- 

0 0 

pond to the cut off frequency of one of the modes of the coaxial 

waveguide, since our formulation is not valid for these frequencies. 

We also note that as w ^ , R(o) ) ^ 0 , thus the TEM mode reflected 

0 0 

field tends to vanish for very high frequency excitations. If 
ka,kb « 1 , it is obvious that for large values of z the higher 
order E^^^ modes are negligible since they are evanescent in this 
frequency range. Since the time behavior of the incident field is 
given, we shall only be concerned with evaluating the integral for 
the reflected TEM mode fields or 



H, (p,2,t) 
<P2 



reflected 

TEM 



1 



M(g)) 









( 5 . 14 ) 



From the values of y^(k) and x^(k) determined in Appendix D, it 
can be adduced that to the order of (ka)^ or (kb)^, [y^(k)]~^ = 0 and 
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[x'"(k)]-2 = -H„ |[j^(ka)] {l-2ika U 



1 1 



n=l l“nl^ 



■)} 



2ka 



(1 - C + Jin 2) 

X e " 



X exp[2 i- K^^^(u) {Jin(1 + — )}du] + 0(ka)^ (5.15) 

x-'n+ /l2 2 

6->-0 V k - u 



By examination of the integral in (5.15) at the extremes of very small 
and very large values of u , it is shown that 

00 

1 K^^^u) {Jln(l + du 

6 7k^- 

is convergent for all values of k . Additionally, the integrable 
singularity at u = k presents no difficulty to evaluation. After 
substitution and a change of variable the integral becomes: 



1 7 ,, x)J„(x) - J,(t x)v^(x)]M,(x)J„(| X) - J,(x)Y„(t X) 

M .. ,b"~ . , x,2 ... xb a ■/ a -.2 I 



0 



[Jo(|x) - Vx)f . [Y„(^x)- Y„(x)]2 



X {jln(1 



1 









-)} dx 



(5.16) 



For small values of x > integral (5.16) can be approximated by 




i 
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I 1 _ 7TX a 7T DX 2 a TT 7TX 2 7T d 

(Itn |)^ 



(1 - 



^x )ln x](^+ 1 ) - ^ )ln ^ 
TT*" **'a ' 7T a 



(in 1)^ 



-) (dx 



This last integral is finite, thus the integral (5.15) is convergent as 
u 0 . For very large values of x , (5.16) is proportional to 

7 ^ + 0(1) + [additional terms]] 

I ^’y[cos(|x- J). cos(x-}) + 0(l)]2+ [sin(^x- J) -sin(x -^) + 0(i)]^} 

X {«,n(l + ^ ^)} dx 



M 



[1 - 



_4_ 

7TX 

- 1 - (2 + 0 ( 1 )) 
7TX ^ 'X'' 



] {«,n(l 



1 






;)} dx 



[1 



7TX 



— (1 + 0 ( 1 )) 

7TX ' 'X'' 



■] (t) dx 



where we have used the expansion of -fn(l + z) with z - A/x ; z « 1; 
z -1 . Additionally the convergence of 



sin mx (^) dx and 



cos mx(l) dx 



M 



has been utilized. We are now able to conclude that the integral 
in (5.15) is convergent as u . 
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Combining all the foregoing, it is apparent that the integral 
appearing in (5.15) is convergent for all values of k and it is appar- 
ent from taking the limit as ka 0 of (5.16) that 

CO 

lim f K*'*(u) (nn(l + — lil — )f du = 0 

Additionally, the term involving the summation in (5.15) is very much 

smaller than one and therefore ignorable. The exponential factors are 

also negligible when large values of z are concerned. Therefore 
^ —2 

[x (k)] in the low frequency limit becomes 

[x'^(k)]'^ - iln ^ [J^(ka)]e"‘''^^ for ka,kb « 1 (5.17) 

a 0 

The reflected TEM mode field is now evaluated by substituting (5.17) 
into (5.14). The results are 

CO 

H. (p.z.t) - JL f ‘iMj (ka) do. (5.18) 

Vo I P 0 

"^reflected TEM -oo 

By the convolution theorem and well known integral (5.8.1) 

f 

(p,z,t) = ■ I ‘1(t'-T) g(x) di ; z£0 (5.19) 

^reflected TEM -«> 



where 
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0 < |t| <7 



9(t) = \ 



M =1 



0 ; otherv/ise ^ 



and 



f = t - - - 

^ c c 



If we again ignore for the moment the frequency restraints and allow a 
unit step function voltage excitation, then the reflected field is 
given by the integral of g(T) . Carrying out the integration produces 
the following: 



H (p,z,t) = 

4>2 

reflected TEM 



^ ^0 Trp iln — ’ 

a 



+ sin 



•1 CT 



T <— \ 



3 ^ ^ a 

— < T < — 

c — c 



I 



J 



T = t - 



c 



£ 

C 



and the electric field is 



Ep (p.z.t) - - H^^(p,z,t) 

reflected TEM reflected TEM 



In the case of a step input the reflected electric field is 
essentially the negative of the incident field for t^^+ ^ 




k 

I 
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(see equation (2.6)). The current reflection coefficient is -J^lka) for 
ka « 1 , which is consistent with transmission line theory. For 
allowable inputs of the form (1 - e”^^^) with T large enough so that 
the low frequency restriction is not violated, the reflected 
fields are given by the above equations, mul tiplied by (1 - e ^) , 
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6. Radiation from the Pulse Excited Cylindrical Monopole Antenna 

The objective of this chapter is to examine the radiation from 
a pulse excited cylindrical monopole antenna in the light of results 
obtained from the preceding study of a coaxial waveguide antenna. We 
have determined that in the far zone the radiation appears to originate 
from the edges of the open ended structure located in the aperture 
plane. The edge starts radiating when the leading edge of the incident 
current pulse on the feed line impinges upon it. Radiation from the 
inner edge is opposite in sign to radiation from the outer edge. 

We begin by considering the approximate far zone fields from 
the annular slot antenna. This antenna has been the subject of several 
authors [28-32]. The structure of the annular slot antenna is that of 
a coaxial line terminated in an infinite ground plane as shown in 
Figure 6-1. Note that the center conductor is not hollow. 

We will follow the approach of Ref. [32] and replace the elec- 
tric field in the aperture by an equivalent magnetic current and 
develop the fields in the half space z >. 0 from this equivalent 
source and image theory. The equivalent magnetic current is 




-a 



z 






Since this source is just above a perfectly conducting ground plane, 
the electric vector potential is given by 
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Figure 6-1. The annular slot antenna 
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2ir b 




ik| r- r' | 



p ' dp ' d4> 



0 a 



where 



e E = - V X F 



0 



and 



|r-r'| = [p^+ p'^ - 2pp' cos(4>-4.') +(z- 



Now It is clear that the electric vector potential and the fields are 
independent of the coordinate 4> , from the cylindrical symmetry of 
the structure. Therefore we can set 4> = 0 without loss of generality. 
We also know that the only nonzero component of F is the <p com- 
ponent due to the symmetry of the equivalent magnetic source. From 
geometry we find 



a , = a cos 4>‘ - a sin <|) 
9 y X 



since 4> = 0 ; a = a . for this formulation. Therefore 

y 9 



F = 




cos 4>' — p'd(J)'dp' 



JkR 

e 



0 a 



with 



R = (p^+ p'^ - 2pp' cos (p' + (z-z')^)^^^ 



From Maxwell's equations and the Lorentz gauge it is an easy 



matter to show that 
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H = [- 7(V • F) * F] 

'‘o 



Since r has only a 4> component which is independent of the angle c|> 



and 



H 



<P 



i (1) F , 
<P 



V • 



0 




p'dp' 



2tt 

I cos ()) ' 
0 



„ikR 

^d<p' 



( 6 . 1 ) 



Equation (6.1) is exact--no approximations have been imposed. If E 
the electric field intensity in the aperture, were known exactly, then 
the integration could be carried out and the result would be exact. 

If we approximate the aperture field as being that of the inci- 
dent field in the aperture, i.e.. 



E'(p.0,o)) = — -- H - 
P p' £n - 



then 



ICJC 






2n 



0 V( 



Hn 



b 2 tt 
a a 0 



ikR 

cos (p' — dp' d(j)' 



( 6 . 2 ) 



In the far zone the distance R becomes 



R = (p^ + p'^ - 2pp' cos 4>‘ + (z- z')^)^^^ 

= (r^+ p'^- 2p'r sin 6 cos 4>' - 2rz'cos 6 + z'^)^^^ 



= r[l + ^ sin e cos p' cos 6 + 



r 
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= r[l - ^ sin 0 cos (|)' ] 
since r » p' and r >> z' , 

gikR 

~R~ "" F 

We have retained more terms of the approximation to R in the phase 
factor than in the amplitude factor, since small variations in R 
cause large changes in the value of the exponent. The amplitude 
of the function is relatively insensitive to variations in R as 
r 00 whereas the phase is not. 

Upon incorporation of the well known expansion [28, p.230] 



* t I * I OO ^ I 

,lkp sine cos J (.if J (kp.si-n e)e'"* 



n=-o 



equation (6.2) becomes 



. . b 2tt ^ 

[ [ cos (})' I (-i)"j (kp'sin 6)e’"‘*’' d(j>' dp' 






ooe 



_o Vja^ e 2 cos^(J)'d4)' J,(kp'sin 6)dp' 



b 27 t 



V(w) ikr f 

—^-5 — -p— J^(kp'sin 9)dp' 

I i 



(joe.V(a)) -ikr 



0 



d(J^(kp'sin 0)) 
k sin 0 



a a 

/S VM e''"" • ''o*''" 

a 



] (6.3) 



sin 0 
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The approximate time behavior of the far zone field follows 



directly from the inverse Fourier transform of (6.3) 



H^(r,e,t) = 



/ 



-io)(t-^) 



V(w)[J^(kb sin 0) 



V ^o (Jln ^)r sin e • 

d ^ 



- JgCka sin 6)]du (6.4) 



When this integration is carried out by the convolution theorem, the 
impulse response of the annular slot antenna is found to be given by 
equation (5.10) and the step response is found to be given by equation 
(5.11). 

From this result it is apparent that in the annular slot antenna 
the radiation also appears to emanate from the edges in the aperture 
plane. The radiation starts when the incident current on the conduct- 
ing surface reaches the edge. One edge is located at p = b and the 
other at p = a . For an incident impulse of voltage the outer edge 
contributes radiated magnetic intensity in the amount 




0 



T ^ sin e 



(6.4.1) 



i 
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and the inner edge contributes radiated magnetic intensity in the 
amount 




T < - — sin 0 
~ c 



- 4sin0< T < -^sin 0 V 
c — — c ( 



T > 4 sin 0 
— c 






with 




(6.4.2) 



We make the assertion at this point that the tip and the base of 
a cylindrical monopole antenna (Figure 6-2) give off radiation in a 
manner identical to that of the annular slot antenna. In this 
case, tiie radiation from the tip must be delayed in time by an amount 
which equals the time it takes light to travel the length of the 
antenna, h/c . We further assume that the currents on the conductors 
are reflected in the same manner as the currents on the conductors of 
the open ended coaxial waveguide. That is, the total current must go 
to zero at the tip and the base. 

In Appendix E we develop the radiated fields of a very long 
hollow cylindrical antenna. Transcribing those results, we have that 
the magnetic intensity of such an antenna is 



H^( r ,0 ,co) 



\u \ sin 0) 

/ 0 V(h)) e 0 ^ 

J y ' - F r sin 0 

^ ^0 sin 0 Jin — 

a 



This result is identical to the frequency domain radiation from the 
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Figure 6-2. The cylindrical monopole antenna 
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tip of the center conductor of the annular slot antenna determined in 
(6.3) above. Since ka « 1 , the particular form of the end does not 
have a significant effect on the field intensities. We therefore take 
these results to be valid for all thin cylindrical antennas regardless of 
the end configuration (with the obvious exception of those wi th end loading.) 

From the above discussion, it is clear that the contribution to 
the radiation from the tip of the center conductor is identical for 
h = 0 and for extrenely large values of h . We assume this result to be 
valid for all cylndrical monopoles with h > 0 . Since the base of the 
antenna is excited at t = 0 , when the incident excitation voltage 
reaches the ground pi ane^ appropriate mathematical factors must be in- 
corporated into the field intensity equations to account for the 
finite time h/c required for the excitation current to reach the 
tip of the antenna, (c is the speed of light). 

Since the electric field at the junction of the feedline and 
the antenna must be continuous, when the flowing charges reflected 
from the tip reach the junction, charges of opposite sign must travel 
back along the antenna [8]. If we consider that the annular thickness 
(b-a) « X , then image theory is applicable and it is readily seen 
that the current relaunched onto the antenna from the base is equal 
in magnitude and opposite in sign to the current incident from the 
antenna. 

No radiation besides the initial radiation due to the current 
incident from the feedline originates from the base of the antenna. 

The current reflected from the antenna tip is cancelled by the image 
current at the z = 0 plane. 
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To determine the total contributions to the far zone radiation, 
use is made of Figure 6-3. Proper accounting must be made of the 
radiation from the tip of the antenna and that reflected from the ground 
plane. It is readily seen from this figure and the above discussion 
that the radiated magnetic intensity in the frequency domain is 



H^(r,0,oj) 



-J 



^0 V(g)) e 



ikr 



*^o (tn sin 9 ° 



[J (kb sin 0) - J-(ka sin 0) 



^ikh(l - cos 0) . ^ikh(l + cos 0) 

} 3 



1 + e 



2ikh 



kh ^ (2n+l) J ; n=0,l,2,3,--- 



( 6 . 6 ) 



The case for kh = (2n+l) will be treated later. 

The time behavior follows from the inverse Fourier transform of 

( 6 . 6 ) 

CO 

H.(r,0,t) = -J- f H.(r.0,o)) e'‘'“^ du) (6.7) 

<P ^ ^ 

-»00 

If the excitation voltage is harmonic time dependent, then 
V(w) = Re[»^ A 6(o)~ w^) ] 



From (6.6) and (6.7), we have 
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Figure 6-3, 



Effect of image plane 
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Hjr.e.t) = 



4> 



^0 {in |-)r sin e 

a 



Re 



w b 



cj h cj h 

— ( 1 - cos 6) i — ^ — ( 1 + cos 6) 



(ij.a 



1 + e 



-iu) (t - — ) 

X {e ° M 



Upon simplification this becomes 



V 



.e.t) =/p 



A cos w^( t- tt) 
0 c 



(Ob 



0 (An — )r sin e 

a 



J«( — — sin e) 



O' c 



(0 a (O h 

sin e)][cos(-|- cos 6)] 



cos 



(o^h 

0 



(o_h 



b f a i (2n+1) | ; n=0,1,2,3. 



( 6 . 8 ) 



This expression shows that the radiation pattern is 



F(e) = 



cob coh coa coh 

J (-^ sin e)cos — [J«(-^ sin e)][cos(-^ cos e)] 

sin e 



(6.9) 



The radiation pattern given by (6.9) can be compared with the pattern 
derived when the antenna is assumed to be very thin with a filamentary 
sinusoidal current distribution [33]. For that case 
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0) h 

cos — — cos (w -cos 0) 

F(6) = . ° (6.10) 

Sin e 

In the limit when the feedline diameters are "thin" (i.e., X » a,b) > 
the ipatterns given by (6.9) and (6.10) are identical. Equation (6.8) 
shows that the radiation along 6 = 0 is always zero, as it must be. 

To get expressions for the cases where kh = (2n+l) ^ , 
n=0,l ,2,3 , . . . , we must consider the fact that ohmic losses are present 
in any real system allowing us to use a complex value of k or we can 
take k to be conplex for the sake of mathematical expediency and let 
its imaginary part vanish at a later stage in the analysis. Thus, 

k = k^ + ik2 ; c 

With complex k the factor 

^ikhd - cos e) . ^Ikh(i t cos e) _ 

j'yTTRli cos kh ^ 

cos[(k.|+ ik2)h cos e] 
coiTFf+”TTc^7^i 

cos(k.jh cos e)cosh(k2h cos 6) - i sin(k.jh cos 8)sinh(k2h cose) 
cos k.|h cosh(k2h) - i sin k.jh sinh(k2h) 

In the limit as k.|h -»• (2n+l) ^ , n=0,l ,2,3, • • • , this expression is 
seen to be finite for k2 ^ 0 and the magnitude is only limited by 
the fact that the conductivity of the material is not infinite 




4 
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lim 

k-| h (2n+ 1 ) 



/ Cos(kh cos e) v 
it' cos kh ’ 
I 



cos [(2n+l cos e] 
-i[sin(2n+l)j]sinh(k2h) 



n=0,l ,2,3, • • • and k 2 = e > 0 



The magnetic intensity becomes 



H^(r,9,t) 




A sin(u t - 

0 c 

(in sin 6 



-0o^“oc^^" e)[cos((2n+l)| cos e)] 
sin[(2n+l)^] sinh(k 2 h) 



The input current to the antenna peaks at the values of kh 
given by kh = (2n+l)^ , n=0,l,2,3,** * , The explanation for this is 
that the monopole can be considered as a simple extension of the feed- 
line. The tip of the antenna would correspond to the open end of a 
transmission line. Any odd number of quarter wave lengths from the 
tip along the line toward the source would correspond to a point where 
the inpedance is a minimum (a point of series resonance) and the cur- 
rent is a maximum. When kh = (2n+l)^ , the feed point is an odd 
number of quarter wavelengths from the open end of the extended trans- 
mission line. 

Taking the inverse Fourier transform of (6.6), we find the 
radiated magnetic intensity for a voltage impulse of excitation 
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1 



H^(r,0,t) = 

impulse 
exci tati on 



/^( 1 

’ % 7r(lln — )r sin 

a 



-) 

e 



sin e) 



; 0 ^ T < ^ sin 6 

2 1 " ^ 



; elsewhere 



T= t-- 




■1 

0 < Q < — sin e 
— n — c 

elsewhere 



f2 ~ t-“ - ^1 - cos 6) 
n c c 



2nh 

c 



i I 



(-1)" 



_ j ./(fsin e)^- 



0 < T < — sin e 
— n — c 

> 

elsewhere j 

T = t -f -^1 ♦ cos e) -2 ?!l (6. 

n c c ' c 



For a unit step of voltage excitation the radiated magnetic 
intensity is the integral of the preceding expression. 
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> T sin 0 
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,^ = t-r-^l.cose) 



(6.13) 



Excitation Voltage 



112 




Figure 6-4. Radiated magnetic intensity at 0 = tt/2 . Ordinate: 



Excitation Voltage 
( vol ts ) 




Figure 6-5. Radiated magnetic intensity at 0 = tt/3 . Ordinate: H = H , ( /iT/e" r £n - sin ^)(vo1ts) 
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t 1 I I L 

2 4 6 



c(t - *"/c) 
b sin 0 



Figure 6-6. Radiated magnetic intensity of annular slot antenna 
Ordinate: H = H ( /u /e r in ^sin 0)( volts) 

(J) q; 0 0 a 

Antenna Dimensions: h = 0 and a = . 8b 
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The magnetic intensity for arbitrary time dependence can be ob- 
tained by convolution of the input excitation voltage with the impulse 
response (6.12) or by convolution of the time derivative of the input 
voltage with the step response (6.13). As expected, the radiated field 
for the infinite monopole is the same as that for the finite length 
monopole for the time interval before the current pulse reaches the tip 
of the antenna. The results, equations (6.12) and (6.13), are easily 
extended to the infinite cylindrical antenna by letting h -> <» . The 
radiated fields for the annular slot are found by letting h 0 in 
expression (6.6). It is also worth observing that the results (6.12) 
and (6.13) are based on the assumption of no ohmic losses as the cur- 
rent pulse travels back and forth between base and tip of the antenna. 
The attenuation can be found using the conventional transmission line 
approach for small attenuations on good conductors. The determination 
of the time domain fields in the case of a lossy conductor is com- 
plicated since the resistance per unit length is proportional to tiie 
square root of the frequency. The incident pulse supplies the total 
of the radiated energy, ohmic losses, reflected energy and energy 
stored in the induction fields. The stored energy of the induction 
fields is the source of energy for the second and subsequent pulses 
launched down the feedline toward the source, the energy radiated after 
time 2h/c, and the ohmic losses which occur after time 2h/c . 

The theoretical results obtained in this chapter are in good 
agreement with the experimental and nunerical work of Schmitt, et al 
[6], Palciauskas and Beam [10], King and Harrison [8], Burrell [34], 
and Lamensdorf [14]. 



-'ll 
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Figure 6-7 is a comparison of the theoretical results obtained in 
this work with the experimental results of Schmitt, et al. [6]. Reference 
[6] states that the second pulse measured in the experiment should be 
moved .74 ns toward the origin and its amplitude should be about 3/2 as 
much as observed. This is necessary to account for experimental devia- 
tions ass oci ate d with locating the receiving probe a finite distance from 
the antenna. In Figure 6-7 these corrections are applied to the second 
and subsequent pulses. 

It is also constructive to include the following observation. 

When a cylindrical dipole antenna is infinitely thin and center fed, it 
is common to consider it as an open ended transmission line that has 
been spread or opened out. A sinusoidal antenna current distribution 
with current nodes at the ends is assuned and the radiated fields are 
obtained on that basis. If it is further assumed that the current at 
the antenna feed point is equal to the current existing at this point 
on an equivalent length of open ended transmission line, then the field 
expressions developed from the assumed sinusoidal current distribution 
are identical to the field expressions obtained in this work in the 
limit of very thin antennas. 

To demonstrate this point, we first observe that the radiated 
magnetic intensity created by a filamentary sinusoidal current distribu- 
tion of ^Qsin k(h - |z|) is [33] 

ile’'*"' 

Lj J_ Q 

“ 2irr sin 6 



[cos kh - cos(kh cos 6)] 
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The driving point current is IpSin kh . The other symbols have the 
usual meanings consistent with their previous use within this chapter. 

It is also common knowledge that the input impedance of an 
open ended lossless transmission line of length h is 



Z = -iZ cot kh 
0 

where Z^ is the characteristic impedance of the line. If we drive 
such a line by a perfect voltage source V(o)), the source current is 
readily obtained by Ohm's law. Since the antenna driving point cur- 
rent is assuned to be equal to the current which would exist at the 
input of an equivalent length of open ended transmission line, we 
equate the driving point current to the source current. Therefore, 



and 



I sin kh = 
0 



-__V(cu) 

-i Z cot kh 
0 



T = 

0 -iZ cos kh 
0 



By substituting this value of I^ into the field intensity equation, 
we derive the general result for very thin antennas, 



H^jj(r,e,w) 



-V(M) 

2frZ r sin 0 
0 



[1 



cos(kh cos 6) 1 
cos kh -* 



( 6 . 14 ) 



Outside of a discrepancy in sign, it is an easy matter to show that 
this general result reduces to equation (6.6) when the limiting form 
for small arguments replaces the Bessel functions of equation (6.6) 
and the characteristic impedance of a lossless coaxial line is 
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substi tuted, i .e. , 




The discrepancy in sign is the result of not specifying a 
voltage reference point and it is not significant relative to the 
results . 

The approximate transient behavior of a very thin dipole 
antenna fed by a transmission line with characteristic impedance 
is found by taking the inverse Fourier transform of equation 
( 6 . 14 ). 
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7. Conclusions 

The results presented in this report provide a comprehensive 
analysis of the open ended coaxial transmission line or coaxial wave- 
guide antenna. Tlie fields radiated from such structures can be com- 
puted without a priori assumption of the aperture field. The current 
reflection coefficient for the fundamental mode is found to be 
-J^(ka) for ka « 1 . 

The conventional approximation for the fields radiated from 
the coaxial aperture antenna gives good agreement with the results ob- 
tained through the Wiener-Hopf method for thin radii in the feedline 
(kb,ka « 1) . 

When the excitation is such that X << b-a, the radiation 
pattern obtained by the conventional analysis method is shown to be 
valid for tiie predominantly forward direction only. We determine that 
the classical radiation pattern obtained by conventional methods should 
be multiplied by the factor |cos 9|'^^ for X 0 . The resulting 
pattern is one in which the radiation at high frequencies is confined 
to the forward direction as expected from geometrical optics predic- 
ti ons . 

A model for the cylindrical monopole or dipole antenna is 
developed based on the observation that radiation appears to emanate 
only from the base and tip of the structure in pulse excitation. This 
model is used to determine the radiated electromagnetic field of 
cylindrical monopole antennas for both harmonic time dependent and 
arbitrary time dependent excitation voltages. 
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The mathematical techniques required to obtain results using 
this model are extremely simple and in most cases hand calculations 
are sufficient. The results obtained show good agreement with experi- 
ment and are useful for both analysis of the radiation field when the 
driving voltage is given and for synthesis of a driving voltage to 
produce a specific radiation field. 

Future research can be aimed at removing the restriction that 
the width of the annular region of the coaxial waveguide antenna is 
very much smaller than the wavelength of the exciting source. The 
results could also be extended to the thick antenna case (i.e., 
kb,ka >. 1). Additional research effort on the transient behavior of 
the cylindrical monopole could be directed toward characterizing the 
effect produced by different end cap configurations, determination of 
frequency limitations, accounting for dispersion and describing the 
effect of higher order mode components of the aperture field on 
the transient radiation. 
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Appendi x A - Analytical Determination that H Modes are not Present 
on Structure 

In this appendix we will show that H modes are not present on 
the coaxial waveguide antenna excited by a TEM mode. 

It is obvious that synmetry precludes excitation of all H modes except 
(n=l ,2,3, ♦ • ♦) modes. 

By the edge conditions [19] (see Figure 2-1) 
we know that the longitudinal current must approach zero. Thus the 
sun of the incident current and the reflected (longitudinal) current 
must be zero. Therefore the reflected current on each conductor has a 
longitudinal component equal and opposite to the incident current. All 
of the reflected current is longitudinal as it must be to completely can 
cel the incident current. To have modes, a circumferential com- 

ponent of current should be present on the structure. Since no current 

source is available to launch a circumferential current, the modes 

on 

cannot be present. 

The above paragraph notwithstanding, let us assume that some 
current is launched in the circumferential direction upon reflection 
of the incident current. It is shown in Chapter 2 that the longitudi- 
nal current in each conductor is equal and opposite. In keeping with 
this, we assume that the reflected currents on inner and outer conduc- 
tors must also be equal in magnitude and directed in opposite circum- 
ferential directions. It will now be shown that = 0 (n=l ,2,3, • • •) 
The applicable scalar Helmholtz equation for this problem is 




I 






;f;5^ 
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9Z 9p p 



(A.l) 



E^(p,z,o)) = 0 at p = a,b ; z _< 0 



E^{p,z,oj) = 0 at p = 0 ^ < 



with edge condition: 



E .(p,z,w) ->-0 at p = a,b ; z ->• 0 
<P 



Fourier transforming (A.l) gives: 



9p ^ P 



(A.2) 



Using the regions of Figure 2-2, with replaced by E^ and 

the asymptotic behavior of E^{p,z,(o) as z -+ ±«> , we find that 
E, (p,a,o)), E (p,a,o)) and E (p,a,o)) are analytic in the strip 



-<t> 



1 



' 4 > 



■4> 



shown in Figure 2-4. Solving equation {A.2) with E^ = E^ , E^ , or 
E respectively, yields 

YO 



= ak^(yp) 



E^^ = bi^(yp) + ck^(yp) 

E = DI^(yp) 



(A. 3) 
(A.4) 
(A.5) 



(A. 3) satisfies the radiation condition and (A.5) goes to zero at 
p = 0 . At p = a , 
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(a,a,o)) 




(a»a,(o) 



E" (a.a.oj) = E" (a,a,cj) = 0 
4>3 4>2 



^pE^ (p,a,o))] 



p=a 3 



p=a 



or 



E^ (a,a,w) = BI^(Ya) + CK^(ya) 



C = 



E^ (a,a,w) - B(I^(Ya)) 



^pE (p,a,o))] 



B + 



From (A. 5) we have 



^pE (p.a.o))] 



= ya 



p=a 



E^ (a.a.w) iQ(Ya) 
I^(Ya) 



where use was made of 



D = 



E. (a>(X}(i)) 
I](Ya) 



Subtracting (A. 11) from (A. 10) and using (A. 8), it follows that 



(A.6) 

(A.7) 

(A.8) 

(A.9) 
(A. 10) 

(A.11) 
(A. 11.1) 
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^pE’ (p,a,o))] - ^pE" (p,a,o))] 

^3 p=a ^2 



p=a 



B + 



Define the left hand side of the above equation, as G~(a) 



E^ (a,a,oj) 

' - iqrfjr * K,ws)i,rYa) 



At p = b : 



E^ (b,a,(o) = E^ (b,a,(o) 

9 ] 92 



E" (b,a,(o) = E” (b,a,u) = 0 

9 1 92 



^pE^^(p,a,03)] 



= (p.ot.o))] 

p=b ^2 



p=b 



or 



E^^(b,a,a)) = BI^(Yb) + CK^Cyb) 

E^ (b.a.oj) = BI^ 
r = ] 



(yb) 



XfW 



■^pE (p,a,o))] 



p=b 



= yb[BI^(Yb) - CK^(Yb)] 



B + 

' K,(yb')- - YbE^^(b,a,«) 



(A. 12) 



(A. 13) 
(A.14) 
(A. 15) 



(A. 16) 



(A.17) 
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(p.a.o))] 



p=b 



+ K (7b) 

-Yb E^^(b,a.o)) 



(A. 18) 



where we have used 



A = 



E^ (b.a.to) 



KjTTa) 



(A. 19) 



Subtracting (A. 17) from (A. 18) and using (A. 15) it follows that 



^pE‘^(p,a.o))] 



)=b V 



p=b 



B 

K^iwy 



Define the left hand side of the above as E~(a) . Then 



E“(a) = 



197 ^ 



(A. 20) 



If we substitute (A. 20) into (A. 12) we get 



G'(a) - E‘(a) 



K^(Yb) E^^(a,a,o)) 
K-|(Ya) ~ l<i (Ya)I-|(Ya) 



(A.21) 



By equating the values of C given in (A. 16) and (A. 9) , it follows that 

+ K,(Ya) ^ Ii(Yb) iT(Ya) 

E (a,a,w) = o iL\ E. (a,a,w) - B[^, ,■ 1 i - u /v^\3 l^i(Ya) 






i^[WT" iqw-* ■"V 



Upon use of equation (A. 20) the above relation reduces to 

+ K,(Ya) + 

E^ (a,a,w) - K (^ - 5 - ) (a.a.w) = E (a)[K^(Ya)I^(Yb) - I^(Ya)K^(Yb)]. 

3 11 



Substitution of this equation into (A.21) gives the companion equation: 



t 




(A.22) 
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G (a) - E (a) 1 ^ 3 ) ~ K,(yb)I,(Ya) 



From the definition of G"(a) we note that it is proportional to 



H (a,a,co) - H (a.a.co), and E“(a) is proportional to H (b,a,o)) 

^3 ^2 

- H (b,a,cj) with the same proportionality factors. By well known 
^2 

boundary conditions we know the surface currents on the outer conduc- 
tors are, respectively: 



parallel plate case is approached. Under these conditions with TEM 
excitation, the currents on the plates can be taken as equal in mag- 
nitude and opposite in direction, with a great degree of reliability. 
Therefore G"(a) = -E“(a) and equation (A. 21) becomes 





In this problem we have assumed that b-a « and therefore the 



G-(a) = 

Ii(Ya)[Ki(Yb) + K,(Ya)] 




(A.23) 



Simultaneously (A.22) becomes 



G-(a) = 

Ki{Yb)[Ii(Yb) + InCYa)] 




(A.24) 



We may factor the denominator of the right hand member in (A.23) and 
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(A.24), using the techniques of Appendix D. 

Ki(Yb)[I-|(Yb) + J-|(Ya)] = N'*’(a)N'(a) 
with N^(a) = N~(-a) as a -+- <» 

with 1 Im a[ > -l <2 

Si mi larly 

Il(Ya)[Ki(Yb) + K^(ya)] = M''’(a) M"(a) 
with M^(a) = M (-a) ''' (la|) as a. 

with I Im a| > -^2 

It is also necessary to point out that G"(a) % as a -► ® . 

This is determinable from the edge condition that H^(p,z,u)) 'v 
as z 0~ at p = a,b . 

Substitution of these factored forms into (A. 24) and (A. 23) 
and algebraic operation yield the entire or integral functions (A. 25) 
and (A. 26) 

(b,a,a>) 

G'(a) M'(a) = (A. 25) 

M'^(a) 

(a,a,o)) 

G'(ci) N'(a) (A. 26) 

N (c.) 

Each side of (A. 25) or (A. 26) is the analytic continuation of the 
other side. Both sides of (A. 25) and (A. 26) are analytic in the strip 
-l <2 < T < l <2 . Since each side of (A. 25) and (A. 26) is bounded for 



-130- 



all a, they must also be constant (Liouville theorem). Examination 

of the asymptotic behavior of the above expressions shows that each 

side is identically zero. (a,a,o)) -»■ 0 as |a| -»■<“ , since 

+ + + 

E (a, 2 ,o)) -> 0 as z -> 0 . Therefore E. (b.a.co) = E. (a,a,u) = 0 . 






1 






Each of the coefficients A,B,C,D can be expressed by a linear combina- 
tion of E^ (b,a,o)) and E^ (a.a.w). This conclusion is evident by 



■<|) 



1 






observation of equations (A. 11.1), (A. 19), and the simultaneous solu- 
tion of 



E^ (a,a,w) = Bl^(ya) + CK^(ya) 
3 

E^^(b,a,o)) = Bl^(yb) + CK^(yb) 



Hence E. , E. and E. are identically zero and E.(p, z,u) does not 
(j)^ ’ (j)2 <J) 

appear in any region of the coaxial antenna structure. 
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Appendix B» Detemi nation of Branch Cuts for y 

In this appendix we will develop the branch cuts in the com- 
plex a plane that define the region of analyticity of Y , a multi- 
valued function of a 

Y = y/a - k 

Our starting point for the work in this appendix will be equations 
(2.9) and (2.9.1); 

p ^ ■ T ■ ° 

dp p ^1 

H. (p,a,oj) -»■ 0 as p -*■ <» (2.9.1) 

4>i 



The solution to equation (2.9) satisfying (2.9.1) was found in Chapter 

2 . 



H 



<^1 



(p,a,d)) = A(a) K^(yp) 



(2.15) 



The asymptotic value of the modified Bessel function [20] is known to 
be 



Ki(yp) 



— e’^*^ as p -*■ «> 



2yp 



Hence 



*1 ( 

H (p,z,o)) « 



A(a) e e da 



as p ->■ 



00 



(B.l) 
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From the Soamerfeld radiation condition we know that the field inten- 
sity must be an outward traveling or evanescent wave for large 



Following the manner of Reference [21, p.20], we note that if 
we requi re that 



may be verified by examination of (B.l). The branch cuts will be 
chosen in a manner which ensures that the conditions (B.2) and (B.3) 
are satisfied. 

The a plane is viewed as a two-sheeted Riemann surface with 
the sheets connected along the branch cut. In each of the sheets the 
function y single valued. The sign of the function on one sheet 
is the negative of the sign on the other sheet. The value of y 
becomes discontinuous only if a branch cut in the a plane is crossed. 

Let us define the branch cut of y such that Re y > 0 in 
the top sheet and Re y < 0 in the bottom sheet. Thus the two sheets 
are joined together by the curve given by Re y = 0 . 

The branch point is obviously k and we will, for computa- 
tional purposes, take k complex 




Re y ^ 0 
Im y £ 0 



(B.2) 

(B.3) 




k = k-| + ik2 , k^ » k2 , k-|,k2 > 0 



We may therefore write 
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= (o + ix)^ - (k^ + ik 2 )^ 

= (a^ - T^) - (k^ - k^) + 2i (ax - k^k 2 ) (B.4) 

Note that Re y > 0 in the entire top sheet, therefore if we 

specify = re^® and y = with 0 _< 0 < Ztt , then 

Re Y > 0 only if |9| < ir . Thus the branch cut must be given by 
2 2 

|e| = TT or Re Y :£ 0 and Im y = 0 . The branch cuts which divide 
the a plane can be easily identified by examination of Figure B-1 
and are seen to be a portion of the curve ox = k-|k 2 . 

Since the branch cut is given by y^ " it follows that 

on the branch cut y= or Re y = 0 . The graph (B-2) depicts 

the branch cuts for the top sheet. 

For the branch cuts in the bottom sheet, it is only necessary 
to use the opposite sign of y everywhere in Figure B-2, 

It is also worthwhile to state that for this branch definition 

•; 0 <}>‘|‘*' *}>2 

Y = /r-|r 2 e where e = — ^ — 

^ 4*] ^ “ "2" » “ ^ ^ 4*2 ^ -^<6<^ 

Y = -i Jk^- or Jk^- o' = iy 



I 



I 



I 

i 

I 






I 



i 
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Figure B-1. Regions of complex a plane 




1 
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Im Y 




Figure B-2. Branch cut for multivalued function y 
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Appendix C. Useful Expressions and Expansions 

In the main body of this work.extensi ve use is made of ordinary 
and modified Bessel functions and Hankel functions. Listed below are 
explicit expressions, asymptotic approximations and limiting forms for 
some of these functions [20,24] of integer order. 



(i) Relations between modified and ordinary Bessel functions: 
- iv iri X Tii 

I^(z) = e J^(z e^ ) ; (-it < arg z < |-) 

K^(z) = j Tii e^ H^^^ze^ ); (-n < arg z < p 



(ii) Wronskians: 









V 



V' 



’v+r 



TTZ 



(iii) Series: 



Jo(z) = 1 






I* 



(j (| z2)3 



( 1 !)^ ( 2 !) 



r 



(3!)‘ 
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1 7 ^ 

2 






Y (z) = f) + C} J (z) + ^ — 5- - y 

o' ' tTI' 2' ^ o' ' TT ^^,)2 2 



{2\y 



(iv) Limiting forms for small arguments (when v is fixed and z -*■ 0) 
J^(z) (^z)V(v+l) ; V f -l,-2.-3,"* 

Yq(z) -iH^^^(z) 'V iH^^^(z) 'x. | Hn z 

Y^(z) -iHj^^(z) iH(,^^(z) Re(v) > 0 

I^(z) % (^zr/r(v+l) ; V f -1,-2, -3,... 



K (z) 'v> -iln z 
0 

K^(z) •>.lr(v)(|)-« 

(v) Asymptotic Expansions (when v is fixed and |z| -> <») 



J^(Z) 


= i|^{cos(z . 


1 

■ I 


c 

1 




(|arg z| < tt) 


Y^(z) 


= ^(s1n(z. 


1 

■ 2 


V7T - — )} ■< 

4 




(|arg z| < tt) 




•. 




(|arg z| 
















K (z) 

V 


„ /n~ ^-z 




(larg z| 


37Tx 
< -Y-) 
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(vi) Derivatives: 

^b\(z)] = -z''K^.,(z) 

(vii) Recurrence Relations: 

I i(z) - I .,(z) = % I (z) 
v-1' ' v+1' ' z v' ' 

K, ,(z) - K ^,(z) = - %. K (z) 
v-r ' v+1' ' z v' ' 

(viii) Other Relations: 

I (z) = I (z) 

-V V 

K (z) = K (z) 

-V V 
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Appendix D. Factorization Procedure 

The purpose of this appendix is to furnish a factorization pro- 
cedure for the Wiener-Hopf kernels which arise throughout this work. 

The factorization procedure used extensively in this work and developed 
in this appendix was originated by Bates and Mittra [23]. The mathe- 
matical development of this procedure is first presented and then it is 
applied to the specific expressions to be factored. 

— Yh 

We begin this appendix with a factorization of e with 

n — ? 

Y = ya - k using the well known theorems of [4]. 

Theorem A : Let f(a) be an analytic function of a = o+ir regular 
in the strip x_ < t < , such that f(a) < C|o| , P > 0 for 

|a| -»• 00 , the inequality holding uniformly for all x in the strip. 
Then f(a) can be decomposed such that 

f(a) = f^(a) + f (a) 

with 

x_ < c < T < x^ 

X < X < d < x_j_ 

where f^(a) is analytic in the portion of the a plane defined by 
X > X and f"(a) is analytic in that portion of the a plane de- 

fined by X < x^ . 



oo+i c 

-00 fi c ^ ~ “ 



dB 



f(a) = 



oo+i d 

— . ^d6 

2iri J B-a 

-oo+i d 



w-r I 

^ f 
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Theorem B ; If Jin g(a) satisfies the conditions of Theorem A, im- 
plying that g(a) is an analytic function of a , which is regular 
and nonzero in the strip t < t < and g(a) 1 uniformly as 
a -*■ ±~ in the strip, then g(a) can be factored such that 



g(a) = g (a) g'(a) 



wi th 



<»+i c 



g'^(a) = exp ^ f dB ; x_ < c < x < 



-co+i c 
o»fi d 



g“(a) = exp - ^ 



~ dB ; x_ < X < d < x^ 



-oo+i d 



The function g^(a) is analytic in that portion of the a plane 
defined by x > x and g”(a) is analytic in that portion of the 
plane defined by x < x^ . The principal branch of the logarithmic 
function is used throughout. It is also worthy to note that if 

and iln g(a) -Jlnjal as |a| ->■«>, Theorem B can still 
The integrals are convergent in the sense 



g(ot) 



1 

■R- 

be applied. 



lim 

T OO 



T+ix 

■ 

-T+i X 



} dB 



If f(a) and g(a) are even functions of a then 



f^(-a) = f'(a) or f^(a) = f'(-a) 
and 
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g^(-ct) = g"(ot) or g"'’(a) = g’(-a) 



The function , which we wish to factor is interpreted as 






^ for or e^^^ 



ri — ? 

~ ^ for . Examine 



«,n g(a) = -yh = -h(a^- k^)(a^- 



The first factor is an entire or integral function of a . Therefore 
we only need to decompose the second factor 1/y . Since y is mul- 
tivalued, we will use the branch cuts for the top sheet of the two 
sheeted Riemann surface of the a plane (see Appendix B). 

From Theorem A, we have 



«>+ic 






d6 



-oo+l c 



^(6-a) 



-k^ < c < T 



(D.l) 



The integration path for equation (D.l) is diagrammed in Figure D-1. 
The integral (D.l) will be evaluated by integration along a contour in 
the upper half plane, which does not cross the branch cut located in 
that half plane (see Figure D-2). By Cauchy's theorem the value of 
the integral (D.l) is 



f*(o) = -Z= 



JT.y. 



ST f 



dB 



7 2rfT j r-K ^ 

P^+ P2 - k‘^(B-a) 



(D.2) 



The contributions to the integral from the semicircular portions of 
the path at 



are zero by Jordan's leirnia. 
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Figure D-1 
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ilmySi 






Re/8 



Figure D-3. 




-144- 



For the limiting case of no losses, k 2 "^0^, the path P-|+ P 2 
is adjacent to the coordinates as shown in Figure D-3. Under these 
ci rcunstances the remaining integral of (D.2) is expressed as the 
five separate integrals which follow; 



1 



^Tii 



{ 



0 



00 



idy 

i \/y^+ (iy-a) 




dx 

(x-a) 



+ lim 



-IT 

f 



i r e^ ® de 



0 

r 



dx 



^ ^ TT 'J 2rke^®+ r^ e^^® [k + re^® - a) ^ -i J^~ ^ (x-a) 



00 

f 



idy 



0 -i (iy-a) 



By reversing the limits on the last two integrals and evaluation of 
the third, we can reduce the above set of integrals to: 



k 

1 f dx 



0 (x-a) x^ 



+ 




dy 

\/y^+ k^ (y+ia) 



A further reduction is achieved by a change of variable. Let 
X = k sin 6 and y = ik sin 6 . Incorporation of this change yields 

ir/2 -i«> 

. i f d6 ^ 1 [ d6 

^ |j ksin6-a ^ ksin6-a 



Note that use has been made of the identity, sin(-iy) = -i sinh y . 
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From Reference [20, p.78] we have the following integral: 



dz 

a+b sin z 



2 a tan b 2 2 

— X — Y'Xr>' 3rctan[ — « — T~YTT^ ’ > b ) 



= \ 



, a tan l + b - (b^-a^)^'^^ 2 ? 

5 5 1 /? 7 5 5 1 /p3 ’ (b ^3 ) 

(b^-a^)*' a tan |- + b+(b^-a^) 



Through this integral we are able to complete the evaluation as fol- 
lows : 



TT 



,u?r?)T7? 



arctan 



r- a tan j + k' 



[ - at< 

TT- 






- arctan 



- a tan |- + k ^ 



(a^ > k^) 



Tr0?^~?r^ 



V 



IT 



An 



-a tan k + (k -a ) 



r-a tan ^ + k - (k^-a^)^"]^ p-a tan k-(k^-a^)^' 

« 2 -?>« I , ■ '1-atan f . k*(k^-c^)‘*. 



(k^ > a^) 



Since tan(-iy) = -i tanhy and lim tanh y = 1 , the above relation 

y CO 

simplifies to: 



( 




\ 



(k^ > a^) 




i 
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The trigonometric identity 



arctan + arctan Z 2 



= arctan[ 



1 + Z^Z2 



•] 



was applied to the top expression. 

Recognizing the relation [25, p.62] 

tan-'z = ^ Zn[ j ; ] 

and applying it, we are able to achieve great simplification 



i 

i 

7TY 






Jln[ 



Y - (a+ k) 

(a + ky 



■] 



(a^ > k^) 
(k^ > a^) 



2 2 2 2 

Note that the same result is found for a > k and k > a . When 
the above is reinserted into (D.2), the decomposition is completed 
and 



f''(a) = 



- [1 + 

Y 



TT ‘‘Y + 



- {a+ k 
(a+ k 



-]] 






TTY (a+ k) - Y"* Y ■ (ot- k)"* 



The factorization of 



g(a) = e 



can now be easily completed. From 



the above we have 
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in g*(a) = * M 



Because g(a) is even in a , it follows that 



and 



iln g (-a) = iln g‘(a) 



ing-w = 



Fi nal ly 



+ lllX 

gV) = g-(-a) = e ' 



(D.3) 



This result is consistent with that of Noble [4, pp. 20-21], and is 
shown to be correct by the following; 






= e 






= e 



-hy 



(D.3.1) 



where In z = i-n|z| + i(}) 

-71 < £ 7T 

The only possible singularity of g^(a) in the upper half plane is a 
branch point singularity. It is obvious that this singularity is not 
present and that g^(a) is regular in the upper half plane, since 
substitution of ±i|y| for y has the same result. Similarly, for 

g'(a) , 
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4> = 2 tan 



-1 liL 

a-k 



Theorem C: [23] Let G(a) be an analytic function of a in the 
region x_ < a < with the following conditions: 



1. G(a) is regular in the region x_ < a < x^ . 

2. G(a) ^ 0 ; G(a) = G(-a) (even). 

3. G(a) 'V' co^ as | Re a| -»■ " where v and h are 

real constants. 

Then G(a) can be represented within the strip by 
G(a) = G'^(a) G~(a) 

where G^(a) is analytic and nonzero in the half plane x > x_ and 
G"(a) is analytic and nonzero in the half plane x < x^ . The 
expressions for G^(a) and G"(a) are given by 



with 



G'^(a) = G'(-a) = (1+ 

{bV) exp[M,n[I^3-l|^]} 



(D.4) 





4 
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(a) 

B^(a) = exp 



1 

Ztti 



00- i d 
-00- i d 



aF(3) 

3(6-a) 



de ; 



|Im a| < |d| ; |t_| > |d| 

(D.4.1) 



(b) 

F(3) = h(3^- - V 2.n(e^- - 2.n c + an G(0) (D.4.2) 

(c) 

2 2 1/2 

Y = (a - k ) ' with branch of Appendix B. 

(d) 



k = k-|+ ik 2 ; k^ » ; k 2 »k^ > 0 and |t_| £ k 2 ; 1 ^2 



We will now show the validity of the above theorem, 
ing that 



B^(-a) = exp 



r_i^ 

2^ 



-id 

f 



-oo_i d 



-aF(B) 
3(6 +a) 




We begin by observ- 



If 3 is replaced by -3 and use is made of the fact that 
F(3) = F(-3) then B^(-a) becomes: 



and 



B^(-a) = exp 



-oofid 

1 f gF(3)d3 

2iTi J 3(3-a) 

oo+i d 



B^(-a) = B'(a) 



We next find the product 
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B^(a) B"(a) = exp 





f qF(3). d3 


Tin J 


6(B-a) J 



c 



where c is the contour shown in Figure D-4. Observe that the asymp- 
totic behavior of the integrand can be used to show that there is no 
contribution to the integral from the vertical segments of the contour 
at |Re b| =o » . By Cauchy's theorem, the value of the integral 

exponent in the above expression is 



2iri I (residues of the integrand, from the poles at 
3 = 0 and 3 = a ) 



Therefore 



-F(0) + F(a) 

hy + ihk + Jin r ^ 

Lg(0)(1-^) J 



B'^(a) B’(a) = 



G(g) exp[hy + ihk] 
^ v/2 

G(0)[1 - 

k^ 



From (D.4), we formulate 



G (a) G (a) = G(0)(1 



2 v/2 . 

B (a) B (a) 
k^ 



exp 





Substituting for B^(a) B (a) and using equation (0.3.1), we find 




m 



i 






■■■* 
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Figure D-4. 



-152- 



6^(a) 6"(a) = 6(a) , which establishes Theorem C. 

Precise knowledge of the asymptotic behavior of G^(a) = 6~(-a) 
as |a| -»■ “ within t > t_ is necessary for successful use of the 
Wiener-Hopf method. As our next task, we will determine this asymptotic 
behavior. Notice that B^(a) is bounded as |a| , therefore by 

letting |a| -»■ <» in (D.4), we find that 



G'^(a) 



= G‘(-a) (ia| ) 



v/2 



lim 

|a|->« 



ih 

7T 



In 



Y+ (a-k) 






oh 1 2a I 

— hr 

7T Ik' 



(D.5) 



for |a| “ within t > t_. In the manner of [23] the factoriza- 

tion formula of Theorem C will now be converted to a form convenient 
for nunerical work. In this development, we concentrate on changing 

+oo-i d 



B^(a) = exp 



1 

2riT 



gF(B) 

B(B-a) 



dB 



-oo-i d 



to a more tractable expression. 

Taking the logarithm and integrating by parts yields 

5,n B'^(a) = I ^B^^ 

-oo_i d 
oe-i d 



= F(s) 

-J-i d 



A(B) tn(^) d6 



where 
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A(8) = 



eh 



6v 



1 dG(e) 



27r(k^-B^)^/^ 2TTi(k^-B^) 27Ti G(B) 



If we let B -*■ -e , then 

-00+1 d 

B*(a) = I A(-6) «ti(2^) 



An 



dB 



oofi d 



This integral will be evaluated on a closed contour in the upper half 
plane by the residue theorem. Therefore 



iln B (a) = - A(-B) Jln(l + |-) dB - 2iTi ^(residues at the poles 

f) ^ of integrand inside 

contour ) 

The path p and the contour are given in Figure D-5. The contribu- 
tion from the semicircular path at |b| is zero. 

We now assume that 

(1) G(a) has at most one branch singularity in the upper half 

2 2 1/2 

plane located at a = +k in the form y = (a - k ) ' . 

(2) G(a) has only simple poles in the upper half of the a 

plane, located at a = +p^^ with | Im p^| > , n=l,2,3,*.* 

(3) G(a) has only simple zeros in the upper half of the a 

plane located at ct = z^^^ with |Im z^^^j > , m-1,2,3,*** . 

Examination of A(b) shows that its only possible poles within 
the contour are those of poles of this function 

occur at the poles and the zeros of G(B) . The residue contribution 



I 
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ilm u 4 






Re u 









Figure D-6. 
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from the poles of G( 6 ) is ( 2 ^) zeros of G(3) it 

is Since G( 3 ) is an even function, when 3 is replaced by 

(- 3 ) the contributions from the poles and zeros become the negative 
of the above. We may now write 



£n 



B'^(a) = - f A(-3) £n{l + |)d3 - I Jln[l+^]+^ <ln[l+f-] 
I ^ n=l Pn m=l ^m 



Set 

u= (kV)l/2^i(32-k2)l/2 

where the branch on the upper Riemann sheet is used (see Appendix B) 




k2)l/2 = 



(k^ 



u2)l/2 



Using this substitution, we have 



£n B (a) = 



A(-3) iln(l + |)d 3 = 



C(u) £n[l 



a 



(k -u ) 



T72-]du 



where 



" h * 2^ ■ 2^ G((k^-u^)^/^)] 



(D. 6 ) 



(see Figure D- 6 ). 

2 2 

From Appendix B we know that on the branch cut (3 - k ) = 0 , 
therefore (3 - k ) is pure imaginary and the variable u is pure 
real. So the introduction of the variable change maps the contour p 
into a contour along the real axis. The contour p' must be in- 
dented around any poles of the integrand which appear on the real 
axis. The inctentations are consistent with pole movement from the 
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real axis if the mediun is assumed slightly lossy. Note also that we 
have accounted for a negative sign by taking p' in the opposite 
direction to the path which p maps into under the variable change. 

Examining C(u), we note that the second term has a pole at 
the origin, and the third term 



1 d_ 
Zui du 



[in G((k^-u^)'/2)] = - 






lG((k2V)V2, 



could possibly have a pole at the origin, and at other locations along 
the real axis, which follows from the fact that u is real. The 
contribution to Hn B^(a) from integrating ov^r the semicircular path 
around the origin is 

^ [Rq" v] ln[l + 

where R is the 
0 

lim u -^in G(k^-u^) (D.7) 
u 0 

1 ^ 

The remaining poles will contribute + I 

n=l 

R (a) = lim (u-p ) ^[Jln G((k^-u^)^^^)] {in[l + — — z— ]} 

where p^ is the nth ordered pole of C(u) on the real axis. 

We may calculate the integral along the straight line portions 



of the path as 



oo 
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I K(u) iln[l + - - - ■ ■ ] du 
6-K)'*' v/k^-u^ 

where 

K(u) = t - 2 ^ [L(U) + L(oe'")] (D.8) 

and 

L(u) = ^[lnG((k^-u^)'/2)] (D.9) 

The integral is interpreted as a principal value type, and it is 
denoted by the bar on the integral sign. The principal value of the 
logarithmic function is understood to have been used throughout. 

In sunmary, we have found that 



tn 



oo 

B*(a) = f K(u) )ln[l + — 2-]du + ln[l + |] 



n=l m=i m 



and the value of G^(a) can be written as 






G^(a) =x/g(0) (1 + n [1 + n [i + ^] 



n=l 



n m= 



n V 



exp 



! f g |} ] - ^ ^ 

X 2 I R„(a)J 



+ — 2 — ]du 

4T7 



(0.9.1) 
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If we perform a similar derivation in which we use B”{a) 
instead of B^(a), as the starting point, ^e result would be; 

oo 

£n B“(a) = i K(u) Jln[l + “ ] du + ^Rq-v) iln[l + ^] 

6^0+ A^-u^ 

-7I KM + I ^n[l - - I ilnCl - 

^ " n=l Pn m=l ^m 

Similar to (0.9.1), it follows that 



R /2 

G-(c) = ,WT (It r) ° n [1 - H (1 - f-) 

n=l Pn nFl ^m 



*7 J Pn<“> 



(D.9.2) 



From (0.9. 1) and (0.9.2), and knowing that G^(a) = G~(a) , we 
recognize that the factorization is: 



R«/2 



a 



G'^(a) = G“(-a) = G(0) (1 + r) ° n (1 + ^) 

^ n=l Pn 



ikh 

T~ 



,;Ot|-)exprM.„[3L^].l 

m=l m L IX/ 

00 

K(u) £n [ 1 + ■ ] du + i I Rpdal )'] 

.A. 2 t ^ ” J 



(0.10) 



I ^ f 



i 
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To conplete the theoretical portion of this appendix, we cite a 
theorem which provides for the expansion of a class of functions into 
infinite products [26, p. 136]. 



Theorem : Let f(a) be a function with its only zeros at a^,a 2 ,a 2 ,*** 

where lim a is infinite and let f(a) be analytic for all a . 

n ^ 00 ” 

Then f(a) can be expressed 



f(a) = f(0) 



f'(0) 

n=l 



{1 



|_,e " 



a ^ 0 
n 



(D.ll) 



If f(a) is both even and entire with its only zeros occurring at 

mr 



±ap±a,,±a„ 



3*‘ ■ ■ 



and lim + a„ = + — , then f(a) can be written 

n 00 



a 

f(a) = f(0) e 




(1 




aa 

mr 



n 

n=l 



(1 + ^) e 
®n 



aa 

mr 



(D.12) 



The exponential factors for the infinite products ensure convergence 
of the separate products. 

Consider 



0 -r> 

n 



za 

mr 



(1 - — ) 

^ mr' 



za 

niT 



^ 1 - 



77 



Thus, 



00 — “22 
n ( 1 - |-) 'Xj 1 - [ as n -► « 

n=l “n n=l nS^ 



1 

Because I converges as n 
n=l n"^ 



“ , the infinite product also con- 



verges as n . 





J 
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We will now apply the foregoing factorization fonmulas to factor 
the expressions: 

{loCya) [Ko(Yb) - K^(ya)] }'^ 
and 

lyvb) [I^(rb) - I^(ya)] )-' 



Taking the top expression first, we define 



’*(<.) = I„(Ya) ■ J„(ay?T7) 

G*'\o) = Ko(Yb) - Kjj(ya) = i A^-o^) 

- H<')(ay?:7)] 

where use has been made of the well known relation between the modified 
and the ordinary Bessel functions [20, p. 370]. 

It is readily seen from (D.12) that 






J(j(ka) 



n 

n=l 



(1 



a 






-)e 



iga 

niT 



n (1 + 

n=l 




iga 

HTT 



where a„ = ±i|a„| ; |a„| k^ ; z„ is the nth ordered zero 

of Jq(z) and z^ -v (n - |-)Tr as n ->■ « . The value of ka must 

be restricted so that la I ^ 0 , therefore take ka « 1 . 

' .n' 

Theorem C is used to factor G^^^(a). From the asymptotic 
value of the modified Bessel functions, we have 
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= o-’/2 e-^l“l /|C(1)’/^ 



as a -»■ ® 



Thus V = - ^ and h = a . From (D.9), 



L<’>(u) 



= ^Un G((/k^V)’/2)> i f[H<’>(ub) -H<’>(ua)]] 



-b Hj^^(ub) + a H^^^(ua) 
H^''(ub) - Hy'(ua) 



From (D.7) 



Rq = ^ ^ L^^^(u) = 0 . 

° u 0 



L(u) has no other poles on 
the real axis, so R,, = 0 
for all n . 



By (D.8) and the analytic continuation H^^^(ze"^^) = H^^^(z) 



K^^^(u) = 



£ 

TT 



^ ^-bH^^^(ub) + aHj^^(ua) 

H^'^(ub) - H^^Cua) 



-bHj^^(ub) + aHj^^(ua) 
H|^^(ub) - H^^^(ua) 



! 



a _ 1 

IT " 2 tT 



2b[Y^(ub)J^(ua) - J^(ub)Y^(ua)] + 2a[Y^(ua)J^(ub) - J^(ua)Y^(ub)] 
[Jq(iA)) - JpCua)]^ + [Y^(ub) - Y^(ua)]^ 



X 
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Factorlzatlon formula (D.IO) gives 






exptif 



, 0 ) 



1/2 



(g- k 



ika 

- -y- 



+ I 
6 ^+ 



£n[l + — ] du ] 



(D.13) 



X(g) = X^(g) X’(g) = {lQ(Ya)[KQ(Yb) - K^Cy^)] } 



-1 



Thus 



X^(g) = 



.X(a) 



0 n=l ^l“nl 



(D.14) 



eX(a) ^5 introduced in the "plus" function and in the "minus" 

function to ensure algebraic behavior of each of these functions as 
|g| -»■«>. Note also that x(~«) ~x(«) since x^(-«) = x'(“) • We 

determine the value of x(ct) asymptotic behavior of the 

denominator of equation (D.14). 

By (0.5) we know the asymptotic behavior of G^^^^(a) is 






(a) 



% 




i 



e 




2a 

T 



as |a| -»■ ® 



in the upper half plane. The asymptotic behavior of the infinite 
product in the denominator of (D.14) can be found from the well known 
property of the gairnia function [4,p.41j, 
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oo r(~ + 1 ) 

„ [, t e-“/“ = 



n-1 



r(|^I* 1) 



Introduce Sterling's formula. 



1 



(X - ^ «_ 1 

r(a) e ^ CL {l + (12a) + 



} as a -*• “ , 

|arg a| < it 



Combining these two relations gives the useful result 

" <. -v: '■(l-fi) | 

a -»c° n=l /cTt 



exp ®1 - C - tn %L]} (D.15) 



In (D.15) C is Euler's constant (0.5772 • . • Applying this re- 
sult we find the asymptotic behavior of the infinite product in (D.14) 
to be 



lim 

Ot ^ 00 



n (1 + 

n=l 









I arg a|<ir 

where Ict^l (n - I-) as n -*■ oo . 

By collecting the above results, we determine that 



x(a) 



{1 

ITT ‘ 



C - tn 



aa 






and that x^(a) = x”(-a) behaves as 
region of the a plane with t > -1<2 . 



+ tn 



2a 



|a| 



1/2 



(D.16) 



as |a| ^ “ in the 
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Tuming to the other expression to be factored, 



K„(Yb)[I„(Yb)- I^(Ya)] 

we define 

f*^’(a) = I^(Yb) - I^(Ya) = J^(b A^-a^) - 
G<^*(a) = KjYb) = i f 



Using (D.12), we have 



= [J^Ckb) 



^ ia/a+bv 

J,0 - TJ^>® 



m=l 






io/a+bx 
iiw' 2 ' 



where a„ = ±i|a„| , 

m ' m' 




and z is the mth ordered zero 
m 

Of J^(z) - J^(|- z) 



The value of kb is restricted to kb « 1 , so that 0 . 

The asymptotic behavior of the zeros of [Jq(z) - Jq(^ z)^ are 
determinable by considering the asymptotic behavior of the Bessel 
function for large argiments and considering that the difference be- 
tween the 'values of a and b does not affect the magnitude. Differ- 
ences between the values of a and b , however, cannot be neglected 
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in phase factors. Through these considerations and trigonometry, we 
find that the most rapidly varying roots are those given by 



m 



(m 



lx 2 TTb 

V a+b 



as m ® ; 



m = 1 ,2 , 3, • • • 



Since 






as a 



then V = - ^ and h = b 



= ^Un(i |)H^l^ub)] = 



-bH^^^(ub) 

H^^^ub) 



= lim u^L^^^(u) = 0 
° u-0 ^ 






njl^ub) H|2)(ub) 

^ 



H^'^Cub) H^'^Cub) 



_ b ri 2 1 

~ ^ TrrrxjT^J 



7T 



J‘(ub) + V‘(ub) 



G<^)*(a) = [|1 H<l)(kb)]'/2 ■ 1 



ikb 

~r 



oo 

I- -f K^^^(u) £n[l + — i^L] dul (D.17) 

6-K)‘*' v/k^- 



Collecting the results, we get 
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V(c.) = Y*(a) Y'(a) = {K „(Yb)[I (vb) - I„(Ya)]}'' 



where 



y^'M = 



g<|)(a) 

. / « ig(a-H)) 

(a) VJ (kb) - J (ka) n (1 + tt^) e 



m=l 






(D.18) 



By (D.5), we know that 



i ob - 1 2a I 

6^^^‘'(a) -V e' as |a| - « 



in the upper half plane and from the asymptotic behavior of infinite 
products (D.15), we know 



lim 

|a|-x» m=l' 
I arg aj <fr 



i / a+b X 
^-T~i 



",0 * Tl^)e“”' - Ic.r'/'' exp{2i|^l -c 



Zn 



g(a^-b) 

2tt 



+ i j]} 



l“ml " FF 



as m ® 



Therefore 



4>(a) 



= [1 - C - tn 



aU+^ 

2 tt 



. iTT . 2b on|2aj -i 



(D.19) 



and Y"^(a) = Y"(-a) -v- as jaj ~ in t > -k^ . 
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Appendix E - Radiation from the Infinite Cylindrical Antenna 

In this appendix we shall determine the radiation from a hollow 
infinite cylindrical antenna. An integral equation for the current on 
such an antenna will be formulated and solved. The fields will then 
be determined from the current. The integral equation considered in 
this appendix was first formulated and solved by Hall^n [35,36]. He 
was interested in the reflected current on a cylindrical antenna, when 
the incident current wave comes from a source that appears to be an 
infinite distance from the end. Although the integral equation is the 
same as that found by Hall^n, the application used in this appendix is 
novel. 

Consider Figure E-1, which is an illustration of the actual 
geometry and its equivalent representation, which is amenable to 
analysis. Since the antenna is made of perfectly conducting material, 
its surface is an equipotential and 

“ "v(t) 6(z) when p=a, 

In the frequency domain this becomes 
£ 2 (^, 0 )) = -v(co) 6(z) 

In view of the rotational symmetry, it can easily be shc 3 wn that 

^ 1 (r) 

E (r,o)) = +io)[A Ir) + V 

r is the vector from the origin to a point on the surface of the 
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Figure E-1. Infinite cylindrical antenna and equivalent repre- 
sentation using gap generator as source 
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antenna and A^(r) is the well known vector magnetic potential, which 
satisfies the relation B = V x A . 

The symmetry of the structure also requires that A^(r) = A^(z), 
thus the axial component of the vector magnetic potential satisfies the 
di f fe ren ti al equati on 

+ k^A (z) = t v((o) 6(2) (E.l) 

We also have from symmetry that A^(z) = A^(-z). 

The solution to (E.l) obtained in a straightforward fashion is 

kvOnle^*^!^! 

Az(z) = C^cos kz + C2Sin kz + ^ 



From the symmetry condition that A(z) = A(-z) and the requirement 
for only outgoing waves on the infinite structure, C-| = C2 = 0 and 



A^(z) = 



2(i) 



(E. 2 ) 



The vector potential will now be determined by an alternative fashion, 
whereby use is made of the free space Green's function and the surface 
current. The general formula is: 



- 





dS' 



with 



|r - r‘ 1= y(z-z')^+ p^+p'^- 2pp'cos(4>-<f)') 



-170- 



In the case of the infinite cylindrical antenna with the obser- 
vation point on the cylindrical surface, this general formula becomes 




-oo 0 



(E.3) 



where 

r = y(z-z')^+ 4a^sin^(^) 



and (f> is taken as zero. Since the vector potential is independent 
of , the particular choice of <}) is inmaterial . Since (E.2) and 
(E.3) are both equivalent expressions for the vector potential, they 
may be set equal to each other. The results are 






2tt .. 

dd)' - ^ 

r p^o) 



ik|z| 



0 



v(„) 



(E.4) 



where -oo < z < oo . 

The equation (E.4) is valid for all values of z , since it is 
based on the mathematical model of the actual antenna structure. 
However, in the actual antenna we observe that the current on the 
feedline after t = 0 , the time when the incident excitation reaches 
the z = 0 plane, is a TEM mode wave. The current wave can be ex- 
pressed as: 

I(z) = 1(0) , 



z < 0 
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1(0) is the magnitude of the current at the plane z = 0 . In view 
of the above, it is justified to represent the antenna feedline system 
with the following set of equations: 




I(z) = 1(0) e"’^^ ; z < 0 (E.5) 



We restrict the analysis to the consideration of TEM mode feedline cur- 
rent only. The function g(z) z < 0 is unknown. The solution to the 

integral equation (E.5) is the current on the infinite cylindrical 
antenna when there is an outward (from the origin) travelling current 
wave on the feedline. 

The current is now expressed in terms of its spatial Fourier 

i n ve rse 



I(z') 



00 




1(a) da 



with 



1(a) = I^(a) + I_(a) 



00 



l'^(a) = f 


I(z‘) 






0 




I"(a) = 1 


I(z') 


~00 





laz ' 



.iaz' 



dz' 



dz' 
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-kpZ kpZ 

Note that |I(z')|£Ae as z-+~ and |I(z')|£Be^ as 
z -*• -«> where k 2 is the imaginary part of the propagation constant. 
1(a) is obviously analytic in the strip -k 2 < x < k 2 of the 
a = 0 + ix plane. I^(a) is analytic in the half plane x > -k 2 and 
I~(a) is analytic in the half plane x < k 2 . 

From (E.5) it follows that 
0 



I-(o) = 



1(0) dz = - i 



m 



(a-k) 



for X < k 2 . Additionally, 

00 

G*(c) = 4- f V(co) dz 



/2tT 



/ 



6 /% 



/2ire^ 



_ , . ""O v(cj) 

■ ' "nr 



for X > -k. 



Upon putting these results back into (E.5), we obtain 



1 

7? 



» 00 2ir . . 

^ dz' [ e"''“^' 1(a) da [ d(})' = 

1 1 0 



-laz 



2tt •’ 

-CO 



G(a) da 



or equivalently 



00 00 



i I I 






il(0) 



-]da 



r” .ifcr 



-00 -CO 



✓77 (a-k) ^ 



d4)' 



= j e-’'“ * G-(a)] do, 



(E.6) 
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We digress at this point to find an equivalent form for 



pikr 

e 



d4>' . It is well known that the free space scalar Green func- 

tion in spherical coordinates satisfies (V^+ k^) G(?,lr') = -6(r-1''), 
and 



G(r.r') = 






In cylindrical coordinates the Green function satisfies 



and 



(v^+ k^) G(p,p',*-*',z-z') = - 6U->‘) «(z-z') 



G(p,p',4.-*',z-z') = -V I 

4tt m=-«o 






Z 2 1/2 

In the above equations y = (ot - k ) ' and p< is the 

1 QC C O ^ 

(greater' ^ * These two Green's functions are equiva- 

lent three-dimensional representations. Therefore, 

j I yvp,)K„(vP,) da 

I r-r I m=-«> j 



and 



2 tt “ 

® d4.' = 2 I^(yp^) K^(yp>) da 



0 



|r-7'| 
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We now evaluate the above equation at p = p' = a and sub- 
stitute it into (E.6). 



00 oo 



1 f 



dz’ f e-'“^'[I*(c.) - 3 do f (ya)K (ra)dB 

J »^(a-k) J ° ° 



-00 -00 



[ [ / 2. ivM+ G'(a)] da 

i ^ ^0 (a+k) 



(E.7) 



Interchanging the order of the integrations on the left hand 
side of the above equation, and recognizing that 

00 

1 dz' = 2„6(o-8) 

-00 

equation (E.7) transforms into: 



OO 00 



2 f da f e-'^^[I^(a) - I (ya)K (ya) 6(a-B) 

J J ^(a-k) ° ° 



dB 



-OO -CO 



= 1 



,->az f 



G'(a)] da 



Integration over the S variable produces the result: 



da {[I^(a) - 



-^r^] Cl.(ra) K (ya)] 
/ST(a-k) ° “ 



r ,- / ’^^O v(oj) . G (a)-| i _ „ 



(E.8) 
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Let ^qCys) factored into the product L^(a) L”(a), 

where L^(a) is analytic in the half plane t > -k, and L"(a) is 
analytic in the half plane t < . 

After substituting the factored form into (E.8) we derive the 
following result: 



oo 

( 

2 



— OO 



da e 



-iaz 



L'(a) {[!*(«) - — ill2L.] L^a) - 

v-^Ca-k) 




G-(a)]) 

(ctk) 



0 



The integrand is regular in the strip -k 2 < x < k 2 of the plane. 
Since it must vanish, it gives the result 



[I+(ot) lUOL] L'"(a) = i )M 

»^(a-k) 2L'(a) ^ (a+k)L'(a) 

2 2 

Multiplication of both sides by a - k yields the integral or entire 
functi on 



L*(o) I*(c.)(a^-k2) - L*(a) = 5^ (^2.^2) 

✓StT 2L (a) 



+ i 




v(to)(g-k) _ 
L‘(a) 



P(cx) 



(E.9) 



The left hand side is analytic in the upper half plane (x > -k 2 ) . 
the right hand side is analytic in the lower half plane (x < k 2 ) and 
each side is the analytic continuation of the other. We divert our 
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attention from (E.9) for the moment to consider L^(a), L'(a) and 
their asymptotic behavior as a -»■ «> in their respective half planes 
of regularity. 

Earlier, we defined the functions 
L''’(a) L‘(a) = iQ(Ya) 



By the methods of Appendix 0, this factorization can be, carried out 
with the results 



L^(a)= L"(-a) = 



n (1 + --i-r) 



laa 

mTT 



m=l 



ml 















ika 

- -Y~ 



oo 

f 

+ - 



6“K) 



K^^^(u) 



J,n[l + 




K^^*(u) = f [1 



2 

TTUa 



1 

jQ(ua) + Y^(ua) 



] 





where 

J„(z) 



z is the mth ordered zero of 
m 



L^(a) 
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as |a| ^ “ in half plane t > -k^ . 



L'(a) 



ot 



- 1/2 



as 



a 



“ in half plane t < l<2 . 



Before returning to (E.9), we note from the edge conditions 
[21] that the current goes to zero at the terminus of the hollow con- 
ductor. As I^(a) and G~(a) represent transforms of current we can 
deduce the asymptotic behavior of these functions for a ~ by 
examining the behavior of the respective currents as we approach the 
ends of the antenna. By considering the antenna length to be finite 
and changing the variable to place the origin at the terminal ends, 
it is easily shown that I^(a) and G~(a) vanish as |a| -*• ® in 
their respective half planes of analyticity. 

Application of the asymptotic behavior of the various func- 
tions to (E.9) shows that 



and 



L^(a) I^(a) (a^-k^) 



i l(0)(g+k) 



a 



as a -»• “ , 



T < -k2 



G' Cq). i£zjil+ i 

L-(a) 2 




v(co) (g-k) 
L'(g) 



< 




as 



g 00 , 



T < k2 



From the extended form of Liouvi lie's theorem, P(g) is a polynomial 
of degree less than or equal to one [4, p.38]. We may therefore write 
P{g) = C^+ C 2 (g-k), where and C 2 are constants to be determined. 

By setting g = k in the left side of equation (E.9), we get 



Addi ti onal ly , 
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i2I(0)k L~^(k) 



X+ Cp = L^{a) l'’(a){a*k) + iUO) L~'(a) , (a-Hk) 

^ (a-k) 

By letting a in the above expression, we find £2 = 0 . Thus 
P(a) is a constant. 

From (E.9) and the value of P(a) just determined, we have 



I+(a) = - _ i2I(0)k L^(k) 

(a-k) »<^(a^-k )L^(a) 

Recall that 1(a) = I^(a) + I”(a), so the transform of the current is 

,(„) = . L*(k) 

^(a -k‘^)L (a) 

. i2I(0)k L*(k) L~(g) 

^ (o.^-k^)I„(ya)K„(Ya) 

Inverse Fourier transformation of the above yields 



I(z) = 



il(0)k L'^(k) 

TT 



^ L~(g) da 

2 (a -k )lQ(Ya)KQ(ya) 



(E.ll) 



It is also clear that if we close the contour in the upper half plane, 

* 1 k z 

the integration produces I(z) = l(0)e ; z < 0 , as required. 
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The vector magnetic potential can now be found from the general 



formula, 



A_ = ^ 



°° 2 .V • I i*> I 

r T/,.\ [ ^ikk-r'l 



_o I I(z') f e 
4tt 27t 



0 



|r-r'| 



d<l)' dz' 



A = 
z 



1p l(0)k L*(k) , f 7 , 

-2—^, 5T J J dz' i 

-0 -OO 



47 T 



ik|r-l^'| 

|r-r'| 



"" I TTT 7 



.-iaz' , » 
e L (g) 



_i (a -k ) JQ(Ya)KQ(Ya) 



da 



iUQl(0)k L‘*’(k) y 

4 tt^ • 

-00 



dz 



' d6[e-’6(2-^')l (ay?37) 



X 



K^(p /a^-k^ )] I 



e L~(g) da 

(a -k )I^(Ya)K^(ya) 



ip^i(o)k L-*-(k) “J 7 



2 tt' 



r 



jdaj 

-00 -oo 



da I (a y?^)K (p 7??) L'(a)e"''^^ 



(a^-k^H^CYa) K^(Ya) 



X 6(a-B) 



iM^I(0)k L‘*'(k) 



e-iaz 

(a^-k^) K^(Ya) 



da 



This equation is valid for p > a , z > 0 . Note that the magnetic 
flux density is given by 
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D - 

4> ■ ■ 9p 

therefore the magnetic intensity is 

il(0)k L'^(k) 7 e'’“^ K^(yp) L"(a) da 

H^(p.2.w) = ^-2 I Y K (ya) (E.12) 

-oo 



Equation (E.12) is exact for an infinite length antenna. If we assume 
that it also holds for an extremely long antenna and employ the 
methods of steepest descent as was done in Chapter 4, the resulting 
equation for the magnetic intensity is given by 



H^(p , z )(i)) 



l(0)e^'^^ 

2rrr 



J^(ka sin 6) 
sin 0 



L"(k) 

L~(k cos e) 



(E.13) 



Earlier we assumed that 1(0) was the incident current which is 
reflected from the discontinuity at the z = 0 plane [8]. The re- 
flected current on the feedline from pulse excitation is knavn through 
experiment to approximate a reflected incident current pulse. From 
the edge conditions for this structure it is clear that the current 
on the outer conductor vanishes at z = 0 . If the current on the 
outer and inner conductors are to be equal in magnitude and opposite 
in direction, then the current on the inner conductor must also vanish 
at z = 0 , the point of discontinuity. We therefore conclude that a 
pulse of current is launched onto the antenna and an oppositely 
traveling pulse is launched back upon the center conductor of the feed 
line. The magnitude of this current is assumed to equal the magnitude 
of the incident pulse. 
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If the incident field is given by equation (2.7) of Chapter 2, 



then 



1(0) = -2tt/-^ 



0 v(m) 



U n D 

0 £n — 

d 



The negative sign is required since the reflected current travels in a 
direction opposite to that of the incident current. Incorporation of 
the above into (E.13) gives the final form of the magnetic intensity. 



..I \ «ikr J^(ka sin 6) i+/u\ 

HJp,z.io) = - 2 ^ -Li^ 

’ V On k Sin 6 



L‘(k cos 6) 



(E.14) 



If 






lim 

ka‘*'0 L'(k cos 6) 



then (E.14) obviously reduces to 



m , 1 = . jM 

V in - ® 

d 



(E.15) 



for an antenna with kh » 1 and ka << 1 . 

To show the validity of this limit we rewrite it from the ex- 
pressions for L^(k) and L~(k cos 6) , 



oo Uii 

n ( 1 + TT§-r)e 
m=l 



Q+(3)(k) 



m=i ''^m' 

ka-^0 n n k cos 6^ ,(-jka cos 9 )/ 6mr) q-h3) / qj g 

' TTF1 

m=l 



'l'(k) 
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exp 



UnD + 



lim 
ka -»-0 



6=0 



7^ 



-]du 



exp 



1 



k2(u) Jin[l + - K-^ - °S -i.. ] du 



L5->0 






lim exp 
ka -»-0 



i 



K^(u) £n[l + J - cos 6 



-] du 



L6-H) 



cos 



6 + /l - (u/k)^ 



This final limit was investigated in Chapter 5 where it was shown to be 
uni ty. 
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